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Abstract 



Using the large acceptance apparatus FOPI, we study pion emission in the reactions 
(energies in A GeV are given in parentheses): 40 Ca+ Ca (0.4, 0.6, 0.8, 1.0, 1.5, 1.93), 
96 Ru+ 96 Ru (0.4, 1.0, 1.5), 96 Zr+ 96 Zr (0.4, 1.0, 1.5), 197 Au+ 197 Au (0.4, 0.6, 0.8, 1.0, 
1.2, 1.5). The observables include longitudinal and transverse rapidity distributions 
and stopping, polar anisotropics, pion multiplicities, transverse momentum spectra, 
ratios (ir + /ir~) of average transverse momenta and of yields, directed flow, elliptic 
flow. The data are compared to earlier data where possible and to transport model 
simulations. 

Key words: heavy ions, pion production, rapidity, stopping, flow, isospin 
PACS: 25.75.-q, 25.75.Dw, 25.75.Ld 



1 Introduction 



The quest to use energetic heavy ion collisions to infer properties of infinite nuclear matter, 
such as the equation of state (EOS) under conditions of density (or pressure) and temper- 
ature significantly different from the ground state conditions has proven to be a difficult 
task, although an impressive amount of data has been obtained using experimental setups 
of increasing sophistication. At incident energies per nucleon on the order or smaller than 
the rest masses, systems consisting of two originally separated nuclei are small in the sense 
that surface effects cannot be neglected even if some kind of transient equilibrium situation 
were achieved. 

In particular, the mean free paths of pions [1], [2], the most abundantly created particles, 
with momenta below 1 GeV/c are neither large nor small compared to typical nuclear 
sizes. Therefore methods of analysis resting on the validity of either of these two extremes 
in the hope to simplify the theoretical description are bound to lead only to qualitative 
success at best. As a result, event simulation codes based on microscopic transport theory 
were soon developed [3] that allowed for multiple elementary collisions as the heavy ion 
reaction proceeds, without however requiring a priori full equilibration. 

An important lesson that was learned in the last two decades was that definite conclusions 
on nuclear matter properties based on a single observable had proven to be premature 
and/or of limited accuracy, aside from not being sufficiently convincing. As an example, 
original hopes [4], [5] to use deficits in pion production relative to expectations based on 
compression-free scenarios were not supported by transport theoretical simulations [6], [7], 
[8]. 

On the other hand transport calculations [9] showed that pion azimuthal correlations, 
'flow', qualify as an observable that could contribute significant constraints on the EOS. 
However, as pion production in the 1A GeV regime (at SIS) is not as copious as it is in 
higher energy regimes and as pion azimuthal correlations turn out to be rather small (an 
effect of a few percent), it was concluded [9] that 'very high statistics and high-precision 
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impact parameter classification are necessary' to exploit the sensitivity of pion flow to the 
EOS. 

This requires the need for large acceptance detection systems capable of registering un- 
der exclusive conditions a much larger number of events than had been possible with 
the Berkeley Streamer Chamber that was operated in the 80 's to obtain pion data at 
the BEVALAC accelerator [10], [11], [12], [4]. The first electronic 47r-detector capable of 
measuring pions was DIOGENE [13] installed at the Saturne synchrotron in Saclay. The 
pion emission data presented in this work were obtained with a large acceptance, high 
granularity device, FOPI [14,15] installed at the SIS accelerator in Darmstadt. Particle 
identification is based on time-of- flight, energy loss and magnetic rigidity measurements 
with use of large volume drift chambers and scintillator arrays. A second large acceptance 
device [16], based on the use of a time projection chamber, was also operated in the nineties 
at the BEVALAC accelerator in Berkeley. In addition, more specialized devices were build 
and used at the SIS accelerator in Darmstadt: KaoS [17], a high resolution magnetic spec- 
trometer and TAPS [18], based on arrays of BaF 2 detectors allowing to identify neutral 
pions and r\ particles via their two-photon decay branches. 

The experimental situation concerning the pion observable before the advent of these 
newer devices was reviewed in ref. [5], some of the more recent particle production data 
obtained in the nineties at the SIS accelerator and also at other accelerators (AGS and 
SPS) covering higher energies have been summarized in ref. [19]. 

Our Collaboration has published pion data before [20], [21], [22], [23], [24], but these were 
more limited in scope and, in particular, did not treat pion flow and isospin dependences. 
Also, we call attention to the fact that the data of ref. [20] (for Au on Au at 1.06A 
GeV) have been revised and should be superseded by the present data (see section 6.1 for 
details). The aim of the present work is to present a more complete systematics of highly 
differential pion emission in heavy ion reactions obtained with the FOPI device, varying 
the incident energy (from 0.4 to 1.9 A GeV), the system's size (from A p + A t = 40 + 40 to 
197+ 197, where A p and A t , are the projectile and target mass number, respectively) and 
the system's isospin, and, of course, the event selection method i.e. the centrality. 

After describing the experimental methods we report on the following observables for 
charged pions of both polarities: 

• longitudinal and transverse rapidity distributions and stopping; 

• polar anisotropies; 

• pion multiplicities; 

• transverse momentum spectra; 

• ratios (7r + //T~) of average transverse momenta and of yields; 

• directed flow; 

• elliptic flow. 

While discussing each of these subjects, we shall refer more explicitly to relevant earlier 
work and compare data where possible. All along we shall also present the results of 
simulations with a transport code showing the degree to which our data can be understood 
on a microscopic level and assessing conditions for equilibration (stopping), sensitivities 
to assumptions on the EOS and the propagation of pions in the medium and searching 
for signals that might give information on the isospin dependence of the EOS. While 
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pions present a probe of hot and compressed matter of high interest in its own right, it 
is also important to have this observable under firm theoretical control as it is a link to 
understanding the production of strangeness under subthreshold conditions where the pion 
emitting baryonic resonances are thought to play an essential role in the collision sequences 
leading to outgoing strange particles, such as kaons. We will end with a summary. 



2 Experimental method 

The experiments were performed at the heavy ion accelerator SIS of GSI/Darmstadt using 
the large acceptance FOPI detector [14,15]. A total of 18 system-energies are analysed for 
this work (energies in A GeV are given in parentheses): 40 Ca+ 40 Ca (0.4, 0.6, 0.8, 1.0, 
1.5, 1.93), 96 Ru+ 96 Ru (0.4, 1.0, 1.5), 96 Zr+ 96 Zr (0.4, 1.0, 1.5), 197 Au+ 197 Au (0.4, 0.6, 
0.8, 1.0, 1.2, 1.5). Particle tracking and energy loss determination are done using two 
drift chambers, the CDC (covering polar angles between 35° and 135°) and the Helitron 
(9° — 26°), both located inside a superconducting solenoid operated at a magnetic field of 
0.6T. A set of scintillator arrays, Plastic Wall (7° -30°), Zero Degree Detector (1.2° - 7°), 
and Barrel (42° — 120°), allow us to measure the time of flight and, below 30°, also the 
energy loss. The velocity resolution below 30° was (0.5 — 1.5)%, the momentum resolution 
in the CDC was (4 — 12)% for momenta of 0.5 to 2 GeV/c, respectively. Use of CDC and 
Helitron allows the identification of pions, as well as good isotope separation for hydrogen 
and helium clusters in a large part of momentum space. Heavier clusters are separated by 
nuclear charge. More features of the experimental method, some of them specific to pions, 
have been described in Ref. [20]. 

2. 1 Pions: reconstructing 4n from FOPI 

The pion data presented in this work are limited to the CDC. In one of two independent 
analyses methods, particle tracking was based on the Hough transform method, HT. Track 
quality cuts were varied systematically and the results extrapolated to zero cuts to achieve 
estimates of the tracking efficiency The relative efficiency of positively and negatively 
charged pions was inferred from studies of the isospin symmetric system 40 Ca+ 40 Ca. An 
alternative method of data analysis with a local tracker, LT, has also been used allowing 
extensive cross checking. This method is documented in ref. [25] and we shall briefly come 
back to it in section 3. 

The measured momentum space distributions of pions do not cover the complete 47r phase 
space and must be complemented by interpolations and extrapolations. In the HT method, 
we filter the data to eliminate regions of distorted measurements (such as edge effects) and 
correct for efficiency where necessary. Since this study is limited to symmetric systems, we 
require reflection symmetry in the center of momentum (c.o.m.). Choosing the c.o.m. as 
reference frame, orienting the z-axis in the beam direction, and ignoring for the moment 
deviations from axial symmetry (see section 9) the two remaining dimensions are charac- 
terized by the longitudinal rapidity y = y z , given by exp(2y) = (1 + (3 Z )/(1 — j3 z ) and the 
transverse (spatial) component t of the four-velocity u, given by u t = fitl- The 3-vector 
P is the velocity in units of the light velocity and 7 = 1/ VI — P 2 - In order to be able to 
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compare longitudinal and transversal degrees of freedom on a common basis, we shall also 
use the transverse rapidity, y x , which is defined by replacing f3 z by f3 x in the expression 
for the longitudinal rapidity. The x-axis is laboratory fixed and hence randomly oriented 
relative to the reaction plane, i.e. we average over deviations from axial symmetry. The 
transverse rapidities y x (or y y ) should not be confused with y t which is defined by replacing 
(3 Z by A = JW+M- 

For thermally equilibrated systems (3 t = \f2f3 x and the local rapidity distributions dN/dy x 
and dN/ dy y (rather than dN/ dy t ) should have the same shape and height than the usual 
longitudinal rapidity distribution dN/dy z , where we will omit the subscript z when no 
confusion is likely. Throughout we use scaled units y = y/y p and u t0 = u t /u p , with 
Up = fip'-fp, the index p referring to the incident projectile in the c.o.m.. In these units the 
initial target-projectile rapidity gap always extends from y = — 1 to y = 1. It is useful 
to recall that in non-viscous one-fluid hydrodynamics many observables scale when the 
system's size and energy, but not the shape or impact parameter, are varied. 

Choosing central collisions of Au on Au at 0.8A GeV as a typical example, we show in the 
upper left panel of Fig. 1 the original tt~ (CDC) data in the (yo-Uto) plane. Proceeding 
from the upper left to the lower right panel of the figure, the next two panels show the 
data after application of a sharp filter and the use of reflection symmetry. 



5 



5 
4 

3 s 3 
2 
1 



Original Data 

_J i I i I i I i L 

-2 -10 1 2 



y 



5 
4 

=r 3 

2 
1 



Symmetry: Reflection 



5 
4 

=r 3 

2 
1 



-2-10 1 2 

y 

~i 1 i 1 i 1 i 1 r 
Fit: Extrapolation 




5 
4 

3 e 3 
2 
1 



"i 1 r 



Sharp Filter 



5 
4 

3 e 3 
2 
1 



J i I i I i I i L_ 

-2-1012 

~i 1 i 1 i 1 i 1 r 
Fit: Interpolation 



-10 12 

y 



f\ ♦ transverse 
* * A longitudinal- 




-2 -10 1 2 

y 



Fig. 1. Distributions dN/dutodyo of pions (tt~ ) emitted in central collisions of Au on Au at 0.8^4 
GeV. The various color tones correspond to cuts differing by factors 1.5. Five panels show the 
evolution of the data treatment (see text), while the sixth panel (lower right) shows deduced 
longitudinal and transverse rapidity distributions. The transverse rapidity distribution with a 
cut on longitudinal rapidity (\uq\ < 0.1), is also shown. 



In the following step each measured phase-space cell dy * du t0 and its local surrounding 



6 



N yu cells, with N yu = {2n y + l){2n u + 1), is least squares fitted using the ansatz 

1 d 2 N 



Uto du to dy 



exp [f(ya,u t0 )] 



where f(x, y) = a 2 x 2 + b 2 y 2 + d\x\y + a\\x\ + b\y + c is a five parameter function. 

This procedure smoothens out statistical errors and allows subsequently a well defined 
iterative extension to gaps in the data. Within errors the smoothened representation of 
the data follows the topology of the original data: typical deviations are 5% i.e. of a 
magnitude that exceeds statistical errors in most cases and is caused by local distortions 
of the apparatus response, thus revealing typical systematic uncertainties. 

The technical parameters of the procedure were chosen to be dy = 0.1, du t o = 0.1, n y = 4, 
n u = 6 (except for the data at 0.4A GeV where n y = 6 and n u — 8). These choices are 
governed by the available statistics and the need to follow the measured topology within 
statistical and systematics errors. Variations of these parameters were investigated and 
found to be uncritical within reasonable limits. 

The smoothened data (middle right panel) are well suited for the final step: the extrapo- 
lation to zero transverse momenta also shown in the figure (bottom left panel). 

The low p t extrapolation procedure was guided by microscopic event simulations (to be 
described in more detail in section 3) that take into account the influence of Coulomb 
effects, which differ for the two kinds of charged pions. 
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Fig. 2. Transverse momentum spectrum of 
negative pions in the reaction Au+Au at 
0.8AGeV. The symbols represent data from a 
simulation with IQMD SM. The solid line is 
a least squares fit (see text) using the data in 
the momentum range (marked 'adjust') from 
0.1 to 0.25 GeV/c. 



Fig. 3. The power parameter used to repro- 
duce the simulated data for positive and neg- 
ative pions as a function of the size (total 
charge) of the system. 



Fig. 2 shows a simulated transverse momentum spectrum for ir~ emitted in central Au on 
Au collisions at 0.8A GeV. It is integrated over all longitudinal rapidities. The range to 
be extrapolated for the experimental data {p t < 0.1 GeV/c) is indicated. The simulated 
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Fig. 4. Measured transverse momentum spectra 
of 7r~ in the reaction Au+Au at 0.8A GeV in 
various scaled rapidity bins of widths 0.1. Cen- 
troids and data scaling factors are given in the 
figure. The data are given by symbols, the solid 
lines represent a smoothened version of the data 
plus extrapolations. 



i i i i 

0.1 0.2 C 

p, (GeV/c) 



0.3 



0.4 



data were fitted in the indicated adjustment range (0.1-0.25 GeV/c) with the function 



the power x and the constant c being two shape parameters, and A" a normalization 
parameter (in the non-relativistic, no-Coulomb regime a 'thermal' fit would mean x = 2.0 
and c = 2mT in terms of a temperature T). The fit function was accepted when it also 
gave a good reproduction of the theoretical data in the low-p f range (even though it was 
not included in the fitting procedure). The best powers x, shown in Fig. 3, turned out to 
be somewhat lower than 2, depending on the system size and the pion charge polarity. 
The dependence on the incident energy was very weak and therefore was ignored. 

The experimental data were then extrapolated using the same (available) adjustment range 
and x-values suggested by the simulation but varying A" and one shape parameter only, c. 
The typical outcome is shown in Fig. 4 for three indicated scaled rapidity intervals. 



2.2 Centrality selection 



Collision centrality selection was obtained by binning distributions of the ratio, ERAT [26], 
of total transverse and longitudinal kinetic energies. In terms of the scaled impact param- 
eter, b = b/b max , we choose the same centralities for all the systems: b < 0.15 (i.e. 2.25% 
in terms of total cross sections), b < 0.25, 0.25 < b < 0.45, 0.45 < b < 0.55. We take 

1/3 1 /3 

bmax = 1.15(Ap + A^ ) fm as effective sharp radius and estimate b from the measured 
differential cross sections for the ERAT distribution using a geometrical sharp-cut ap- 
proximation. ERAT selections do not imply a priori a chemical bias. Autocorrelations in 
high transverse momentum population, that are caused by the selection of high ERAT 
values, are avoided by not including identified pions in the selection criterion. 



Nexp(-pf/c) 
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3 Simulation using the IQMD transport code. 

In the present work we are making extensive use of the code IQMD [27] which is based on 
Quantum Molecular Dynamics [28]. Many interesting predictions and revealing interpre- 
tations of the mechanisms involving pions have been made [29]- [34], [9] with this model 
before our data were available making it almost mandatory to confront it now with the 
present measured data. 

One of the motivations for numerical simulations of heavy ion reactions with transport 
codes is the possibility to investigate the effect of the underlying equation of state, EOS, 
on the experimental observables without relying on the restrictive assumption of local (and 
a fortiori global) equilibrium followed by 'sudden' freeze-out of all elementary hadronic 
collisions. In between collisions, nucleons are propagating in mean fields, the nuclear parts 
of which correspond in the limit of infinite matter to well defined zero temperature EOS. 
These EOS can be chosen to be 'stiff' or 'soft', as characterized in terms of incompress- 
ibilities K = 380 MeV, respectively 200 MeV, where K = 9p 2 d 2 (E /A) /dp 2 near p = p , 
the saturation density. We shall henceforth label these two options, HM, respectively SM. 
The M in HM and SM stands for the momentum dependence of the nucleon-nucleon inter- 
action. IQMD incorporates a phenomenological Ansatz fitted to experimental data on the 
real part of the nucleon-nucleus optical potential. In the IQMD code pions are produced 
by the decay of the 1232 MeV A baryon resonance and may be reabsorbed exclusively 
by forming a A again. The dominance of the lowest nucleonic excitation in the 1A GeV 
energy regime has been directly demonstrated in proton-pion correlation studies [35] , [36] , 
[37]/ 

While in between collisions, the A baryons are propagating in the same nuclear and 
Coulomb fields as the nucleons, pions are feeling only the Coulomb potential. For sim- 
plicity, and in order to limit the number of input variables we 'shut off' the poorly known 
(isospin) symmetry potential in these exploratory calculations. Thus 'isospin effects', if 
any, would have to result either from the isospin dependent NN cross sections (imple- 
mented in the code), e.g. a 'cascade' effect, and/or the Coulomb fields. As in the model 
pion production and absorption is strongly connected with A baryon production and ab- 
sorption, the physics of A propagation in the medium becomes important both for the 
observed final number of pions and the flow pion 'daughters' inherit from their 'parents'. 
If not otherwise stated we have used in our calculations the scheme of ref. [8] for the A 
baryon mass distribution and width. 

Some recent implementations of transport theoretical codes for heavy ion collisions in the 
1A GeV regime contain more advanced features than those just described. Without at- 
tempting to be complete, we mention the works of Bao An Li and coworkers [38] and of the 
Catania group [39] which implement various isospin dependences of the mean field. Other 
codes have made progress implementing non-equilibrium aspects of the local densities [40] , 
better accounting for off-shell effects when particle (resonance) lifetimes are short and/or 
the collision rate is high [41] , reassessing the in- medium cross sections for the NN — > NA 
reaction [42], etc. . The code of Danielewicz [43] includes nucleonic clusters up to mass 
three and has been used in our context here to clarify the role of expansion on the pion 
observables. Very recently, the relativistic code UrQmd has also been used [44] to study 
SIS energy reactions, although its original design [45] was directed towards much higher 
energies. Most of these codes include the influence of many higher baryonic resonances. 
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Fig. 5. Simulated reduced impact parameter distributions for Au+Au collisions at 1.5 A GeV 
using the global observable ERAT for event selection. The two peaks correspond to nominal 
centralities bo < 0.15 and 0.45 < bo < 0.55, respectively, as indicated by the vertical lines and 
the double arrows. Histograms (crosses) correspond to unfiltered (filtered) data. 

One of the advantages of a QMD type code over BUU (Boltzmann-Uehling-Uhlenbeck) 
implementations might be its capability to induce finite (nucleon) number fluctuations 
without resorting to additional assumptions or recipes. For the simulation of the experi- 
mental situation it is of paramount importance that there is an accurate 'centraiity match- 
ing', i.e. the centraiity criterion of the experiment (here the observable ERAT) and its 
realistic fluctuations should be simulated. We vary b uniformly between and 1 and add 
the events (about 50000 per system-energy) weighed by bo- Fig. 5 shows b distributions 
resulting from ERAT selection in two of our standard centraiity intervals. The nominal 
bo intervals are also shown, as well as the effect of applying a filter that takes into account 
the geometrical and the threshold limits of the apparatus. As can be seen this filter does 
not have a dramatic influence on the distributions, and hence its details are not critical. 
However, the effect of using ERAT, rather than the experimentally elusive bo, is not neg- 
ligible, even with a perfect detector. For the energy range of interest here, we find that 
the alternative selection method, binning charged particle multiplicities, yields less sharp 
bo distributions as long as bo < 0.40, but is competitive for more peripheral collisions. 
In particular, multiplicity binning avoids the high bo tails beyond bo = 0.6 visible in the 
figure. Multiplicity selected data are available, but will not be shown here. ERAT has 
the advantage that it does not require the simulation to reproduce the degree of nucleonic 
clusterization, a difficult (still) task for transport codes (see ref. [46]). Finite number ef- 
fects deteriorate the b 'resolution' when lighter systems, such as Ca+Ca, are studied, an 
additional reason to try to be realistic when simulating the experiment. 

We have also used the IQMD code with a more technical aim: we have introduced it [25] 
as an event generator for a GEANT based [47] Monte Carlo simulation of our apparatus 
response to better assess tracking efficiencies and losses due to geometrical limitations. 
This alternative independent data analysis [25] using a local tracking method, LT, instead 
of the Hough-transform based tracker, HT, was found to yield 47T pion multiplicities in 
good agreement (10% or better) with the method outlined above (e.g. see Figs. 17 and 20 
in section 6). 
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4 Rapidity distributions and stopping 



Before looking at special projections in momentum space it is useful to display the full (2- 
dimensional) distributions dN / dutody$. One such distribution was already shown in Fig. 1 
for Au+Au at 0.8A GeV. To see the evolution of the topology with incident energy - i.e. 
from OAA to 1.5A GeV - two more such distributions are shown in the upper panels of 
Fig. 6. 
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Fig. 6. dN/dutodyo distributions for pions (ir~) emitted in 6y < 0.15 collisions of Au+Au at 
0.4^4 GeV (left) and 1.5A GeV (right). Top panels: data, bottom panels: simulation. The grey 
(color) tones differ by factors 1.5. 

The pion (here ir~) sources are centered at midrapidity with no readily evident memory 
of the initial conditions (i.e. around \yo\ = 1). Comparing the contours at OAA GeV 
with those at 1.5A GeV, one notices that in the scaled units chosen for these plots, the 
distributions are significantly wider in both dimensions at the lower energy. The data at 
0.4^4 GeV span more than twice the full rapidity gap. At least two properties involved in 
pion production are not expected to scale: a) the nucleonic Fermi motion which should be 
relatively more important at the lower energy and b) the decay kinematics of the parent 
A (and other) baryonic resonances. 



The simulation with IQMD (see the lower panels) qualitatively reproduces the higher 
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'compactness' at 1.5A GeV, but is less expanded in the longitudinal direction, especially 
at the lower energy. 




Fig. 7. Various excitation functions for bo < 0.15 collisions of Au+Au. The symbols represent 
from top to bottom the stopping observable vartl, the variance of the longitudinal rapidity, the 
variance of the transverse rapidity. The corresponding observables from the IQMD simulation are 
represented by solid (dashed) lines for SM (HM). The dotted curve is the measured [48] stopping 
for nucleons. 

A more quantitative assessment of these features can be made in terms of longitudinal, 
dN/dy 0z , and transverse, dN/dy 0x , rapidity distributions deduced from the reconstructed 
47r data. Due to apparatus limitations, and in order to keep the technical definition of 
the scaled variances of these distributions strictly constant over the full range of system- 
energies and centralities, we define the variances o~ 2 (yo z ), resp. <7 2 (y 0x ), of these distri- 
butions in the finite interval |yo| < 1-8 for pions. Further we define the ratio vartl = 
(Vox) / & 2 (yoz) which we shall loosely call 'degree of stopping' or just 'stopping'. In a non- 
relativistic purely 'thermal' interpretation this ratio would represent the ratio of trans- 
verse to longitudinal 'temperatures', which, if different from one, would then imply non- 
equilibrium. In Fig. 7 we show measured excitation functions for <r 2 (?/ 0z ), <7 2 (y y0x ) and of 
vartl (symbols indicated in the figure). 

The scaled variances are seen to decrease significantly and steadily with incident energy, 
as expected from the qualitative discussion of Fig. 6. The 'pionic stopping' is consistently 
below one, but decreases less rapidly with energy. We also include from our earlier work [48] 
the 'nucleonic stopping' (dotted) for comparison: it decreases faster with energy than its 
pionic counterpart. Some of the difference between pions and nucleons results from the 
fact that the observed pionic momenta result from a convolution of the excited nucleon 
momenta with the decay kinematics, a convolution that is expected to lead to a more 
homogeneous (isotropic) final momentum distribution. A more subtle effect would be that 
different hadrons witness on the average different collision histories in a non-equilibrium 
situation. 

These effects should be taken care of in a microscopic simulation. The solid lines in Fig. 7 
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are the prediction of IQMD SM which follow the measured transverse variances amazingly 
well, but show a flatter trend for the longitudinal variances coming closer to the data at the 
highest energy and underestimating the measured values at the low energy end, confirming 
again trends seen qualitatively in Fig. 6. As a consequence the 'stopping' is also overesti- 
mated especially at the low energy end. The mean field (SM versus HM) has a modest, 
but not negligible, influence: the stiffer EOS leads to a general increase of the stopping 
(uppermost dashed line) by roughly 5% caused by the correlated decrease/increase of the 
longitudinal/transverse variance. 
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Fig. 8. Pion stopping, vartl, as function of the number of participant nucleons, A par t for collisions 
at 1.5^4, 1.0 A and 0AA GeV (top to bottom). Left panels: data, right panels: IQMD simulation. 
The two branches correspond to two modes of varying A par t: 1) central (bo < 0.15) collisions 
varying the total system size, marked 'vary A sys \ and 2) collisions of Au+Au varying centrality, 
marked 'vary bo'. 

Naively, if stopping is incomplete, the degree of stopping is expected to depend on the 
system's size characterized by the total number of nucleons, A sys . Our data comprise the 
systems Ca+Ca (A sys = 80), Ru+Ru (A sys = 192) and Au+Au (A sys = 394) where we 
roughly double the nucleonic size from one system to the next. In order to compare, we keep 
the centrality in scaled impact parameters constant (bo < 0.15). In terms of 'participants' 
one can also vary the (participant) system size A part , by changing the impact parameter bo, 
keeping A sys constant. (In order not to interrupt the flow of arguments here, we refer the 
reader to section 6 for a definition and discussion of A part .) With the concept of participant 
size, we are then in a position to establish with our data a summary of the pion stopping 
systematics, based on the two methods just discussed, see Fig. 8 and its caption. In the 
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figure we have plotted averages for both polarities, 7r + and ir , as no significant difference 
was observed. 

The main message of this systematics (which is reproduced semiquantitatively by the 
simulation) is that there is no unique dependence on A part , but also a dependence on 
either the shape of the participant volume (which depends on the collision geometry, 
hence 60), or on the presence of more or less 'spectator' matter in the early stages of 
the reaction, an interpretation we favour in the present energy regime. The idea is that 
pions, being fast moving light hadrons (notice in Fig. 6 some pions are seen to freeze out 
with 4-5 times the incident nucleon velocities u p ) are able to penetrate into the spectator 
matter before the latter has left the neighbourhood of the participant zone. These pions 
are then rescattered experiencing a partial longitudinal reacceleration. This results in a 
smaller apparent stopping, an effect that increases with the size of the spectator matter, 
as suggested by Fig. 7. These pion (or A baryon) rescatterings in spectator matter also 
influence the anisotropics of azimuthal emissions [29,30] to be discussed later. 
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5 Anisotropies 



In thermal model analyses of pion data (see for example [49]) it is often assumed that 
pion emission is isotropic in the c.o.m. in order to extend to 4-7T the measured mid-rapidity 
data. However, already early studies [50,51,2] of inclusive reactions have reported devia- 
tions from isotropy: usually the polar angle distributions have a minimum near 90° c.o.m.. 
More recently, the TAPS collaboration has observed [2,52,53] anisotropic it emission in 
asymmetric heavy ion systems at subthreshold energies and interpreted the inclusive data 
in terms of a 'primordial' pion emission followed by final state interactions (rescattering 
and absorption). The averaging over impact parameters and the presence of asymmetrically 
distributed spectator matter in the chosen reactions complicate the microscopic interpre- 
tation of such data. In the present work, the fact that the stopping observable vartl (see 
previous section) was found to be less than one, indicates that isotropy is not fulfilled even 
in the most exclusive central collisions, despite minimal amounts of 'shadowing' spectator 
matter. 

To quantify anisotropy it is useful to pass from the coordinates (y z ,p t ,<f)) to spherical 
coordinates (p,0,(f>), where 9 is the polar and the azimuthal angle. This coordinate 
change does not introduce new information, of course, but allows to define observables that 
show deviations from isotropy more directly and more sensitively, especially when vartl 
is only slightly below 1. The distributions dN/duQ t dy such as those shown in Fig. 6 were 
transformed to spherical coordinates duly taking into account the Jacobian. Fig. 9 shows 
the projection onto the polar angle axis for Au+Au at 1A GeV for various centralities. 

We define as 'anisotropy factor' Af the ratio 

+i 

Af = J f(x)dx/2a 
-l 

where 

/(cos 9) = a + a 2 cos 2 9 + a 4 cos 9 

is least squares fitted to the data. Af is the factor with which a measurement restricted 
to 90° (c.o.m.) or to mid-rapidity, has to be multiplied to correct the measured yields 
for deviations from isotropy. For the cases shown in Fig. 9 we find Af = 1.20 ± 0.02 (for 
b < 0.15) and A f = 1.44 ± 0.04 (for b = 0.5 ± 0.05). 

In terms of our Af the anisotropy measured [54] for central collisions 40 Ar + KC1 at 1.8 A 
GeV , Af = 1.19, compares well with our values of 1.21 ± 0.03 and 1.22 ± 0.03 for 6% 
central collisions of 40 Ca + 40 Ca at 1.5 and 1.93A GeV, respectively. 

The polar angle anisotropies cannot be explained in the framework of thermal models. 
The simulation with IQMD, however, does a reasonable job for this observable, aside 
from the absolute normalization to be discussed later (see right panel of Fig. 9). High 
polar angle anisotropies suggest that a significant fraction of the observed pions must 
be produced and emitted close to the system's surface where they can escape without 
rescattering. One could conjecture that some, potentially more thermalized, pions (or 
nucleonic resonances) created deeper inside the system did not reach the detectors because 
of in-medium absorption. 
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Fig. 9. Polar (c.o.m.) angle distributions of ir~ mesons in the reaction Au+Au at 1A GeV. 
The solid lines are least squares fits of the three-parameter function ao + ai cos 2 9 + 04 cos 6. 
Left panel: Experimental data. Open squares (blue): 60 < 0.15, triangles (red): 60 = 0.50 ± 0.05. 
The data have been normalized to 1 at 90° to emphasize the differences in shape. Right panel: 
Comparison with IQMD for 0.25 < b < .45. The ratio IQMD/FOPI is also shown (Thick solid 
curve) . 
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Fig. 10. Excitation function of the anisotropy factor for Au+Au collisions with 60 < 0.25. The 
(black) crosses are the measured data (with systematic errors), (blue) diamonds and (red) squares 
are predictions from IQMD SM and HM, respectively (with statistical errors). 

More systematic details reveal differences between the data and the simulations. While 
stopping (Fig. 7) shows a weakly decreasing trend with increasing energy, the measured 
excitation function of Af looks rather flat (within the indicated systematic errors), see 
Fig. 10, while the IQMD simulations always underestimate Af and predict a rise with 
beam energy which is slightly more pronounced with the softer EOS. 

The observables of anisotropy, Af, and of stopping, vartl, are not equivalent since they 
quantify different aspects of momentum space population. However, they are related to 
the degree that high transparency that would yield vartl <C 1 should be characterized 
by Af ^> 1. The two branches seen in the stopping when varying either the size of the 
system or the geometry, Fig. 8, are also seen with the anisotropy observable, Fig. 11. 
IQMD simulations reproduce these features qualitatively. A somewhat puzzling feature is 
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Fig. 11. Anisotropy (averaged over tt~ and ir + ) at 0.4, 1.0 and 1.5^4 GeV incident energy as 
a function of A part . Full symbols and solid lines represent the data, open symbols and dashed 
lines represent the simulation IQMD SM. The lines are linear least squares fits added to guide 
the eye. As in Fig. 8, there are two branches, one varying A par t via centrality (bo) binning in 
Au+Au collisions, the other via changes of the system size (A sys ), keeping the centrality constant 
(bo < 0.15). 

the fact that the measured data seem to suggest that the anisotropy grows slightly with 
A part when the centrality is kept high (b < 0.15). 

In Fig. 12 we take a more differential look at anisotropy displaying for bo < 0.25 collisions 
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Fig. 12. Anisotropy of it in Au+Au collisions (bo < 0.25) at 1.5, 1.0 and 0.6^4 GeV as a 
function of the scaled momentum no- Full (black) diamonds: data, open (blue) diamonds: IQMD 
SM, open (red) squares: IQMD HM. 
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the dependence on the scaled momentum w = (p/m 7r )o {m^ pion mass) for three different 
incident energies. The simulations with IQMD SM and HM are also shown. Again, IQMD 
underestimates Af systematically, shows a weak mean- field sensitivity and some tendency 
to saturate around uq > 3. Except maybe at 1.5^4 GeV, the data do not show such 
saturation: Af rises approximately linearly within the uq range shown and converges to zero 
at zero momentum. This is expected on general grounds and indicates some consistency 
in our extrapolation to low p t which affects primarily Af(uo) values below Uq = 1. 

For inclusive data it appears that the momentum (or kinetic energy E^ in ) dependence 
of Af shows a maximum near uq = 2.8 (or E^—lbO MeV at E/A = 0.8A GeV) as 
demonstrated in Fig. 13 which reproduces the early data of Nagamiya et al. [51] for 
Ar+KCl at 0.8A GeV. It is interesting to note that the authors concluded (in 1981) 'these 
features cannot be explained by any conventional theoretical model! . In the figure we have 
included a calculation with IQMD for a similar system-energy, Ca+Ca at 1A GeV, which 
is seen to reproduce fairly well the observed features. Inclusive data are generally more 
difficult to interprete in a definite way. Due to trigger biases used to enhance central 
collisions and minimize background reactions [20] the present data cannot be directly 
compared with inclusive data such as those shown in Fig. 13. 
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at 90° for inclusive reactions. Circles: Ar+KCl 
at 0.8^4 GeV [51]. Squares joined by a smooth 
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6 Pion production 



6.1 Experimental trends 

In this section we present our pion multiplicity data and compare the results with those 
of earlier work, as well as with those of transport model calculations. 

In the earlier literature some authors reported on pion production in terms of pion mul- 
tiplicity per 'participant', M(tt) / 'A part . The motivation for this choice came from the ob- 
servation, probably first reported in ref. [10], that in relativistic heavy ion reactions this 
reduced multiplicity did not appear to depend on the size of the system if the incident 
energy was kept fixed. Like impact parameters, 'participants' (or 'spectators') are not di- 
rect observables, they are not defined rigorously and are estimated using recipes that may 
vary with the authors. This limits the level of accuracy with which such 'reduced' data 
from different experiments can be compared. In experiments that measure multiplicities 
eventwise there are two ways to determine the number of participants: first, one defines all 
nucleons to be participants that have momenta outside the Fermi spheres around target 
and projectile momenta, and, second, one calculates from the 'impact parameter' the geo- 
metrical size of a straight-trajectory, sharp-geometry overlap zone [55] where the 'impact 
parameter' is estimated from some global event observable (ERAT, charged particle mul- 
tiplicity, etc) as described earlier. In the present work, using the second method and the 
observable ERAT, we estimate a (90 ±5%) 'participation' for our most central, bo < 0.15, 
sample. 




Fig. 14. Dependence of the participant frac- 
tion on the scaled impact parameter bo in the 
sharp-cut geometrical model. 



b 

Although the calculation of A part with the sharp-cut geometrical model [55] is a simple 
mathematical exercise, we show the result in Fig. 14 to raise the awareness of an approx- 
imate fact: the participant fraction A part /A sys decreases linearly with the reduced impact 
parameter bo if one confines oneself to bo < 0.8. Within this limit, it also does not matter 
much if the geometrical model allows for a diffuse [20], rather than a sharp surface. 

For inclusive measurements, i.e. pion multiplicity data taken without a 'centrality bias', 
Apart = l/4v4 sys is often assumed. 

To our knowledge no single experimental apparatus measures all three isospin components 
of the pion. Assuming that Coulomb effects on pion production are very small, isospin 
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symmetry implies that M(vr + ) + M(tt°) + M(7r - ) = 3 * M(tt°) = 1.5 * [M{it + ) + Af(7r")]. 
If only M(7T~) has been measured, assumptions must be made to obtain the full pion 
multiplicity, M(tt) = f * M(tt~). Using for the following only the data for the most 
central collisions, our excitation functions for the total reduced pion multiplicity are simply 
excitation functions of the quantity [M(ir~) + M(n + )] * 1.5/ (0.9A sys ). 




beam energy (A GeV) beam energy (A GeV) 

Fig. 15. Excitation functions of the reduced pion multiplicity. The left panel compares the 
present data for Ca+Ca with the data (rescaled by a factor 1.1) of Harris et al. [12] for Ar+KCl. 
The right panel compares the present data for Au+Au with the La+La data [4]. 

The two-panel figure shows our measured pion-multiplicity excitation functions for the 
systems Ca+Ca (left panel) and Au+Au. They are compared with Streamer Chamber 
data for Ar+KCl [11,56,12] and for La+La [4], respectively. The smooth curve (a second 
order polynomial) has been fitted to the FOPI data only. The agreement between the two 
sets seems to be excellent. 

At this point we have to remark that the present new Au on Au data are in conflict with 
our earlier publication [20] concerning pion production at 1.06A GeV. A reassessment of 
our older data has lead to the following conclusion: due to a too low setting of the potential 
voltage of the CDC the Chamber was not fully efficient for low ionizing particles (ref. [20] 
had reported the very first application of the CDC). 

A few further comments on the comparison with the Streamer Chamber data are useful. 
In the La+La experiment [4] a 384 scintillator hodoscope covering angles 9i a b < 18° was 
added to the Streamer Chamber and the projectile-spectator angle window was defined 
'to be the region centered about beam-velocity Z=l,2 fragments containing 90% of the 
charged particles in minimum bias'. An efficiency factor, estimated from a cascade code 
simulation, was applied. The authors also mention a correction (approximately 14%) to the 
pion multiplicities to account for track losses and misidentifications. Finally, the negative 
pion multiplicities were multiplied by a factor f=2.35 to account for ir + and 7r° emission. 
From our 7r~/7r + systematics we know that the factor depends on the incident energy (we 
estimate f=2.06 at .4A GeV and f=2.43 at 1.5A GeV for La on La). 

Concerning Ar+KCl, the number of 'participating protons' and the ir~ multiplicities are 
conveniently given in a Table [11]. Fragments with projectile velocity in a 4° forward 
cone and positive tracks (around target rapidity) with laboratory momenta pi a t < 200 
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MeV/c were counted as 'spectators'. All 'participant' tracks were assumed to be singly 
charged. This is a good approximation at the highest energy. From the geometrical model 
we estimate a participant total charge of 28 (for a 180mb trigger) in reasonable agreement 
with the tabulated values [11] for beam energies above 1A GeV. The system Ar+KCl 
(A sys = 77.25, Z sys = 36 and hence N sys /Z sys = 1.146) is not a strictly isospin-symmetric 
system. Using the simple-minded 'isobar' formula [5] one obtains M{n~) / M{ji + ) = 1.255, 
i.e. a number significantly different from 1 and implying f=2.70. Again, using our energy 
dependent systematics, we get factors f=2.55 at 0.4A GeV and f=2.74 at 1.5A GeV, rather 
than f=3 mentioned in [12]. From the tabulated 'proton participants' and M(n~) in [11] we 
get M (ir) /Ap ar t values very close to those shown in [12] if we apply the isobar model f factor 
to the n~ multiplicities and the factor A sys /Z sys to the proton participant multiplicities. 
In our figure we show that we get perfect agreement with our data if we boost the Harris 
data by 10%. No correction for the Streamer Chamber response is mentioned in the three 
publications [11,56,12] which appear to pertain to the same experiment. 

At the GSI/SIS accelerator site pion data have also been obtained using two other experi- 
mental setups. The pion multiplicity data of the TAPS collaboration have been summarized 
in [49]. Some (unpublished) data (for E/A=lGeV) from the thesis work of Wagner [57] 
measured with the KaoS setup, are listed in a Table of ref. [19]. 

A common feature of these experiments is that they are triggering on the particle(s) of 
interest (charged pions in the KaoS case and two-photon candidates for a 7r° decay in the 
TAPS case), within an angular range that is substantially limited if compared to the large 
acceptance devices (but can be varied by moving the apparatus). As a consequence they 
are not measuring pion multiplicities event by event, but rather differential cross sections 
in a restricted part of momentum space, preferably 'around mid-rapidity', under more or 
less exclusive conditions. To convert the, apparatus-response corrected, data to reduced 
47r pion multiplicities one must divide by a nuclear reaction cross section, a r which in 
reference [49] was taken to be a r = irrl(A^ 3 + A^ 3 ) 2 with r = 1.14fm. Further, an 
extrapolation to phase-space outside mid-rapidity has to be made and possible biases in 
the triggering modes have to be corrected for if inclusive data are wanted. One asset of 
the TAPS method is that momenta are measured all the way down to zero. In Fig. 16 we 
compare the present Ca+Ca data and some of our earlier data [22], [23] with the inclusive 
TAPS data for neutral pions (multiplied by f=3) [49] for the systems Ar+Ca, Ca+Ca, 
recalling that the FOPI data are high centrality data. For completeness we also mention 
here TAPS data taken at much lower energies (Ar+Ca at 0.18A GeV [58]) well outside 
the range of the figure. 

The published mid-rapidity TAPS data have been obtained under the assumption of 
isotropic emission. Using our anisotropy information (previous section) we have boosted 
these data by a factor 1.25. If this is done, a satisfactory agreement with the FOPI data 
(and the old Streamer Chamber data) is obtained, the TAPS data having somewhat more 
straggling around the smooth trend inferred from our data (solid curve). 

The figure also shows inclusive data for C+C, [59], again, corrected for anisotropy, and 
for central collisions of Ni+Ni [22,23]. The two dashed curves represent the Ca+Ca curve 
(solid) rescaled by factors 0.82 and 1.58. Clearly, there is some system size dependence 
of the reduced pion multiplicity. We note that the C+C system corresponds to A part = 6 
(inclusive data), the Ar+Ca (TAPS) data to A part = 20 and the FOPI data to A part = 72 
and 100 for Ca+Ca and Ni+Ni, respectively. 
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Fig. 16. Excitation functions for the reduced pion multiplicities in reactions of Ca+Ca (smooth 
solid line identical to the one shown in Fig. 15), Ar+Ca and Ca+Ca [49], Ni+Ni [22], and C+C 
[49]. The respective symbols are indicated in the figure, for further comments see text. 
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Fig. 17. Multiplicities of ir + and ir~ mesons in 96 Ru + 96 Ru reactions at 0.4^4 (left) and 1.5^4 
GeV versus the reduced impact parameter &o- The data points marked with larger full symbols 
are from an earlier publication [24]. 

Participant size dependences can also be shown by varying the collision geometry keeping 
the system fixed, as demonstrated in Fig. 17. The straight lines are linear least squares fits 
demonstrating (within error bars) the approximate linear dependence on the (reduced) 
impact parameter, which implies, as shown in Fig. 14, an approximate linearity also in 
terms of A part . This figure also contains data from an earlier publication [24] of our Col- 
laboration which used the alternate analysis method ('local tracking') briefly described in 
subsection 3. 

To obtain a quantitative evaluation of the A part dependence we confine ourselves first to 
the present data, thus avoiding systematic differences between various experiments and 
hopefully profiting from the smaller point-to-point errors. In the upper left panel of Fig. 18 
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Fig. 18. Reduced pion multiplicities and power law fits. The data points marked pp (or p) and 
CC (or C) are from refs. [60] and [59], respectively. The data in the left lower panel are TAPS 
data [59,49,61]. KaoS data [57,19] are shown in the right lower panel. See text for further details. 

we show the participant-size dependence of M n / 'A part at an incident beam energy of 1A 
GeV. 

It is of some interest to try to join up such data all the way down to the nucleon-nucleon 
system A sys = 2. This raises the question [19] how to define the number of 'partici- 
pants' in such a reaction. Specifically for the pp reaction, we follow ref. [19] and use 
(ai n /a to t)[M(7T + ) + M(tt")] * 1.5/2 where the multiplicities are obtained from the compila- 
tion (Table 4) of [60] and the ratio of inelastic to total cross section is set at 0.4 (for a 1A 
GeV incident beam). The resulting data point is marked pp in the figure and plotted at 
Ap ar t — 2, simply because 2 is the minimum number of 'participants' for any reaction. The 
second 'non-FOPI' point, marked CC, has been determined by the TAPS collaboration [59] 
for inclusive reactions (A part = 6) and represents 1.25 * 3 * M(ir°), the factor 1.25 account- 
ing conservatively for the (unknown) polar angle anisotropy, as discussed before. All other 
data are from the present work for various centralities plus a Ni+Ni point after [21] (for 
Apart = 90, Fig.l). In the Figure the solid line is a power law fit, M(it)/ 'A part = c* A par t T ~ 1 , 
using only the FOPI data (i.e. A part > 20). We find r = 0.86 ± 0.04, significantly different 
from one. (If it was one, the plotted reduced multiplicity would be constant.) This fit, 
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somewhat surprisingly, is nicely compatible with both the pp and the CC data. Including 
all the data in the fit lowers the uncertainty of the r parameter by a factor 2. 

The power law curve is shown again in the other three panels. In the upper right panel of 
Fig. 18 we show only the data points for the most central collisions (except for the pp and 
CC data). The power law fit made to the more extensive set of data in the upper left panel 
is also optimal for this more limited set of data. The advantage of a limitation to very 
central collisions is that one is less dependent on the definition of A part since A part ~ A sys . 
At the higher energy, l.bA GeV, the same type of data, also shown in the figure, are again 
fully compatible with an A T dependence, although r = 0.90 ±0.02 is somewhat closer to a 
pure 'volume' dependence. In deriving power law descriptions of system-size dependences 
one hopes to distinguish between 'volume' and 'surface' effects in a way reminiscent of mass 
formulae, although the relative accuracy of mass measurements is multifold better (and 
therefore the results more convincing). Fitting the data with c v A + c s A 2 ^ 3 instead of cA T , 
i.e. again with a two-parameter Ansatz, one can describe the 1A GeV data equally well as 
shown in the figure (dashed curve). One finds c v = 0.096 ± .011 and c s = 0.186 ± 0.033. 
With this description one can tentatively interprete the size dependence in terms of a 
decreasing surface to volume ratio: the fitted coefficients suggest a 52% surface effect for 
C+C decreasing to 21% for Au+Au (both for b < 0.15 collisions). 

In principle, one could also add an isospin term, again following mass formulae. The only 
unambiguous information in our data comes from a comparison of 96 Ru + 96 Ru with 96 Zr 
+ 96 Zr. However, for the total multiplicity of pions (adding both charges) we found no 
significant difference within errors. 

From the smoothened trends described by the power law fits we can deduce that the 
reduced pion yields in the reaction Au+Au around 1A GeV are lower by a factor of about 
0.85 with respect to the Ni+Ni system. In Fig. 2 and Table 3 of ref. [21] we had concluded 
that this factor was about 0.53. As mentioned earlier, the partially inadequate operation 
of the CDC in the experiment leading to the results published in ref. [20] forces us to 
retract this number and most of the quantitative aspects of the conclusions in this early 
publication. 

The TAPS data for 1A GeV [61,59,49] are shown in the lower left panel together with the 
pp point. We have applied a factor 1.25 to these mid-rapidity data to correct them for the 
polar anisotropies discussed earlier. The data from the earlier publication [61], (red) full 
diamonds, for the systems Ar+Ca, Kr+Zr and Au+Au, have been later revised [49], black 
open triangles, to account for a trigger 'bias towards centrality'. However, the Ar+Ca 
point was not changed. These data were obtained in coincidence with the Forward Wall of 
FOPI. The C+C point was obtained in a different setup not involving the FOPI apparatus. 
As can be seen from the figure, after the revisal, the three TAPS points (triangles) appear 
to be aligned, but are significantly below the solid power law line (the dashed curve is our 
power law curve down-scaled by a factor 1.58). Note also that if the correction is due to a 
centrality bias the revisal also leads to a shift along the A part axis. Trying to join up the 
three revised points to the CC and pp data, one has to introduce two kinks in the overall 
TAPS curve. 

Pion multiplicity data from the KaoS Collaboration are cited in [19] and taken from an 
unpublished part of the thesis work of Wagner [57]. The two data points for 1A GeV 
incident beams, cited as being inclusive, for Ni+Ni, resp. Au+Au, are shown in the lower 
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left panel as taken from Table 2 of ref. [19]. Also shown in the right lower panel of Fig. 18 
is a reassessment of these data points: first, a reaction cross section in line with our 
geometrical scaling (r = 1.15) and the TAPS procedure (r = 1.14) and, second, the 
application, again, of the anisotropy factor Af = 1.25. When this is done, the Ni+Ni point 
completely, and the Au+Au point marginally, agree with our power law fit. 

It is of some interest to mention here that 7r-nucleus reaction data from the meson facto- 
ries [1] have also been described in terms of power laws. In particular, 'true' absorption [62] 
is characterized by r = 0.75 ± 0.05 for (200-400) MeV/c pion momenta. However, the ini- 
tial states of the 'inverse' (absorption) reaction cannot be compared simply with the final 
states in the heavy ion induced pion emission. Indeed, in the latter case, a 'bulk' or volume 
term for a strongly interacting particle like the pion (or its parent A baryon) can only be 
understood if one invokes a collective expansion mechanism allowing pion absorption to 
stop or 'freeze out' globally at some time in the evolution of the reaction. This mechanism, 
presumably, is absent in 7r-nucleus reactions. 
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Fig. 19. Excitation function for mid-rapidity 
7r + mesons in central Au+Au collisions. The 
three data points at the higher energy are from 
ref. [63]. The data are for a centrality of 5% (in 
terms of the reaction cross section) within a 
rapidity interval |yo| < 0.25. 
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The data of our present work for central Au+Au reactions join up smoothly to AGS 
data [63] if taken under the same exclusive conditions. This is shown in Fig. 19 for 7r + 
mid-rapidity data. Included in the figure is a least squares fitted second order polynomial, 
in terms of \n(E/A), that reproduces the six data points with an average accuracy of 6%, 
close to the point-to-point systematic errors of both collaborations. We conclude that there 
is no significant 'kink' in this excitation function, as one might have expected to observe 
if the gradual passage to 'resonance' matter (i.e. with a significant degree of nucleonic 
excitations) were to be associated with a fast increase of entropy. 



6.2 Comparison with transport models 



The old Ar+KCl and La+La data obtained at the BEVALAC accelerator with a Streamer 
Chamber [12,4] seemed to indicate within error bars that the multiplicity per participant 
was independent on the size of the fireball at a fixed incident energy. The authors [4] con- 
cluded that 'pion production is a bulk nuclear-matter probe rather than a surface probe' 
and 'is unaffected by the expansion phase'. As a consequence they deduced by comparison 
with thermal model expectations that the nuclear EOS ('missing energy') 'is extremely 
stiff. Probably the first publication putting into question the sensitivity of pion multiplic- 
ities to the EOS was ref. [6] in which a transport code based on the Boltzmann equation 
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Fig. 20. Left panel: excitation function for total pion multiplicity in central Au+Au collisions. 
For the data at 0.4, 1.0 and 1.5^4 GeV two different trackers have been used, HT and LT. The 
measured data are compared to a simulation using IQMD SM (full squares). Right panel: ratio 
theory/experiment (HT). Full squares: IQMD SM with the systematic uncertainties of the data; 
open triangles: IQMD HM. 

and including the mean field and Pauli blocking effects was used. Later, Kruse et al. [7] 
using a different code confirmed that they could not determine the EOS from the data 
of Harris et al. [12,4], since they found only a 10% effect when changing from a stiff to a 
soft EOS. In 1986 Kitazoe et al. [8] reported that they were actually able to reproduce the 
Harris data 'without introducing the nuclear compression effect. 

Fig. 20 shows a comparison of 47r integrated pion multiplicities (1.5 * (tt + + with 
the predictions of IQMD for central (b^ < 0.15) Au on Au collisions. The ratio theory 
to experiment is found to be 1.22 ± 0.08 for the soft EOS (SM), a number that holds 
independently of the tracker (LT or HT) used. The HM (stiff) version predicts at the 
highest energy (1.5^4 GeV) a drop of pion multiplicities by about 10% and even less at 
lower energies. Danielewicz [43] showed that the 'missing energy' conjectured to be the 
compression energy was actually taken up by the collective (radial) flow energy generated 
during the expansion which has a high degree, but not perfect, adiabaticity (and undergoes 
the associated cooling and memory loss). 

In ref. [64] eight different transport codes were compared: the predictions for pion mul- 
tiplicities in reactions of Au on Au at 1A GeV differed by a factor 1.6 if the 'standard' 
versions were used. Some of this disturbing finding could be resolved by 'unifying' the treat- 
ment of the A baryon lifetime in the medium which was found to strongly influence pion 
production. Performing an IQMD simulation with a 'phase shift prescription' [65,66,64] 
instead of the prescription proposed in ref. [8], we found that the pion multiplicities for 
central Au+Au collisions were reduced by (26 — 30)% to values somewhat lower than the 
experimental values. As we shall see later, the phase shift prescription also influences the 
predicted pion flow, although not quite as significantly. 

In ref. [42] a realistic estimation of in-medium modifications was performed. Perhaps the 
most spectacular and interesting effect originates from a significant modification of the 
'elementary' free cross sections for the NN — > NA reaction caused by the drop of the 
baryonic Dirac masses in the medium predicted presently by a number of theoretical 
models (see for example [67]). Within the quasiparticle picture this affects the kinematical 
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Fig. 21. Ratio of theoretical and experimental pion multiplicities for central Au on Au collisions 
versus beam energy. Crosses and full circles are calculations from ref. [42] with and without 
in-medium corrected cross sections, respectively. The open squares are IQMD calculations with 
free cross sections. The shaded band reflects the uncertainty of the experimental data. 

phase space factors in front of the square of the transition matrix element determining 
the cross sections. Although, due to partial thermalisation and hence a setting in of the 
back reaction via detailed balance, this does not translate linearly to the observed pion 
multiplicities. The effects predicted [42] are still significantly larger than the experimental 
uncertainties. 

Medium effects are of special interest in this context since A baryon resonances decay in the 
medium due to their very short lifetime and undergo a cycle of several regenerations [33] . 
At present it seems that the theoretical situation regarding these matters is not settled. 
This is illustrated in Fig. 21 where we plot the ratio of calculated pion multiplicities to 
experimental data for Au on Au at various incident energies. The ratios were evaluated 
using BUU results from ref. [42] for the so-called NL2-set2 combination of parameters 
(see the original paper for details) which seemed to come closest to a preliminary version 
of our data. Compared to the experimental uncertainty band indicated in the figure the 
difference between the calculation fully accounting for in-medium modifications and the 
'standard' calculation with free cross sections is significant at the two higher beam energies. 
For comparison we also plotted the 'standard' (in the present work) IQMD predictions. 
The standard versions of IQMD and BUU make different predictions. At 1.5 A GeV the 
standard IQMD value is close to the in-medium corrected value of the BUU code. 



7 Transverse momentum spectra 



The main characteristics of pion emission transverse to the beam axis have been described 
in section 4 in terms of the variance of the transverse rapidity distribution which was 
compared with the variance of the longitudinal distribution. In this section we briefly 
compare transverse momentum spectra in more detail with, first, data from the KaoS 
Collaboration [68] and, second, with the output of IQMD calculations. 
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7.1 Comparison with KaoS 



To convince ourselves that the filtered part of our data (upper right panel in Fig. 1) that 
we choose to analyse is not affected by p t dependent distortions in the p t range that we 
are interested in (i.e. excluding both very low p t which we cannot measure, as well as 
very high p t which are rare), we have made a comparison with KaoS data [68] for the 
reaction Au+Au at 1.5 A GeV which is shown in Fig. 22. Notice the comparison is shown 
on a linear ordinate scale. To avoid the influence of assumptions and possible problems of 
aligning the centralities, we have chosen a relatively large sample (the 1200mb innermost 
centrality) at fixed laboratory angles indicated in the figure. The agreement is good as 
far as the shape of the spectra is concerned (the KaoS data have been rescaled to agree 
with the FOPI data on an integral basis). Below 0.1 GeV/c the FOPI data are affected 
by apparatus response. 

7.2 Comparison with IQMD 

Attempts to describe transverse momentum spectra in terms of thermal model parameters 
are generally successful, see for instance [22,23,24]. Here we shall not go through the 
thermal exercise, but are rather interested how well microscopic calculations that are 
free of ad hoc adjustable parameters are reproducing our data. In Fig. 7 we have shown 
that the variances of the transverse rapidities were rather well reproduced in Au+Au 
collisions. We normalize the calculations to the experimental data, since we have already 
discussed (see Fig. 20) the integrated pion multiplicities and are now interested in the 
shapes of the spectra. The outcome of the comparison is shown in Fig. 23 for the reaction 
Au+Au at 1.5 A GeV and b < 0.15. The, somewhat different, shapes of both the 7r + and 
7r~ mesons are well reproduced in the p t range shown. The level of accuracy on which 
the comparison is meaningful in the figure is determined by the systematic errors in the 
data and corresponds approximately to the size of the symbols. We recall that IQMD 
takes into account the effect of the Coulomb fields on the pion trajectories, but does not 
include (yet) higher resonances which might influence the spectra at higher momenta. 
With this partial success in hand, we can go back to the details of the collision history 
that a microscopic code is able to furnish. Rather than being a spectrum resulting from 
the decay of resonances moving locally around with a single temperature, the model tells 
us [31] that the measured spectrum is a superposition of pions from various A baryon 
generations representing different stages of the collision with spectra varying in 'hardness' 
as time proceeds. This is a key to memorize some effect of the early compression stage 
in contrast to a completely equilibrated freeze out scenario and is probably the reason 
why an, albeit small, effect of the mean field on both the yields and the stopping can be 
seen. As we shall find out later, the mean field effect also influences the small azimuthal 
anisotropics. 
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Fig. 22. Laboratory momentum spectra, da / {dpia^d^l) , of emitted pions at 40° (left) and 60° 
(right) for Au+Au at GeV. The KaoS and the FOPI (triangles) spectra are compared for 
the same centrality selection. 




Fig. 23. Transverse momentum spectra of ir~ and ir + mesons in central collisions (bo < 0.15) of 
Au+Au at 1. 5 A GeV. The symbols represent IQMD simulations, the solid lines are generated from 
the smoothened representation of the measured data. The part to the left of the line at pt = 0.1 
GeV/c is extrapolated. The spectra are integrated over the longitudinal rapidity \jjq\ < 1.8. 
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8 Pion isospin dependences 

The 7T~ and the 7r + mesons are members of an isospin triplet. A comparison of observables 
connected with differences between the 7r~ and the 7r + mesons offers therefore, in principle, 
the possibility to explore isospin effects on the reaction dynamics and perhaps also on the 
more fundamental issue, the isospin dependence of the EOS. 

In this section we show our results for the ratio of average transverse momenta and the 
ratio of yields. Another important observable, isospin differences in the pion flow, will be 
treated separately in section 9. 

8.1 Average transverse momenta 

In the simulation with IQMD, pions are assumed to propagate in Coulomb fields. As these 
change sign with the charge of the meson, higher average transverse momenta, < p t >, 
are predicted for 7r + relative to tt~ mesons. This is shown in Fig. 24 for central collisions 
of the systems 40 Ca+ 40 Ca, 96 Ru+ 96 Ru, and Au+Au as a function of beam energy. The 
calculations predict an increase of the ratio < p t > + / < pt >~ with the total charge of 
the system, as one would expect, and an approximately linear decrease with the beam 
energy. The experimental data follow these trends, confirming, for this observable, the 
dominant influence of Coulomb fields, presumably after the strong interactions have frozen 
out. In principle this effect can be used to infer the size of the freeze-out volume. For an 
application of this idea using analytical formulae we refer to [69,24]. Alternatively, since 
the microscopic simulation describes the transverse momentum data well, one could derive 
the freeze out conditions from the detailed intermediate output of the code. Here we rather 
invert the argument and conclude that the code is realistic as far as freeze out densities 
are concerned. 

A closer look into Fig. 24 suggests that the decrease with beam energy is somewhat faster 
at the lower energy end and has a tendency to saturate at the higher energies. This effect 
is marginal however, in view of the experimental uncertainties, which are primarily caused 
by the necessity to extrapolate the measured data to zero p t . 

8.2 Ratio it-/it + 

The usefulness of the 7i^/n + yield ratio for investigating the isospin dependence of the 
EOS has been advocated recently [70,71]. The lessons learned in connection with total 
pion multiplicities (see section 6) leave room only for a cautious optimism. 

It is argued [70] that both extremes of modelling, first chance collisions dominated by the 
A-resonance mechanism (the tt~ /tt + ratio is then w (N/Z) 2 [5]), as well as equilibrium 
statistical models, via the difference p n — p p in the neutron and proton chemical potentials, 
all depend on the isospin asymmetry densities quantified by p n j p p or p n — p p (p n , p p are the 
neutron, resp. proton densities). Under the influence of the heavy ion dynamics changes 
in local asymmetries may be induced that in turn create different n~ /n + generations. In 
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Fig. 24. Ratio of the average transverse momenta of ir + and ir~ mesons. Shown are the excitation 
functions for the systems Ru+Ru (top panel), Au+Au and Ca+Ca (lower right). The IQMD SM 
predictions (full symbols) are joined by linear least squares fits to guide the eye, the open symbols 
represent values inferred from our 4tt extrapolated data. 

particular, if the EOS is stiff against isospin asymmetry at high densities, the local energy 
density will tend to minimize by pushing the neutrons out of high density areas (in neutron 
rich systems) leading to a local lowering of the 7r~ production mechanism. 



In the various panels of Fig. 25 we show a summary of our measured tt~ /tt + ratios and 
compare with IQMD predictions. Briefly, one observes a decrease of the /n + ratio with 
incident energy (upper left panel) which is qualitatively also predicted by IQMD. However, 
in this panel, and more so in the other three panels, it is seen that while IQMD is doing 
a perfect job at l.hA GeV, also when (N/Z) is varied, it clearly underestimates the pion 
ratio for large (N/Z) at OAA GeV. The right lower panel repeats the comparison at OAA 
GeV, but for the filtered data, leading to the same conclusion and thus showing that the 
extrapolations to 4tt are not responsible for the discrepancy. 

The linear (N/Z) dependence at the higher energy instead of the expected (N/Z) 2 depen- 
dence of the A-resonance model (the model used in IQMD) can be partially understood 
when realizing that the copious pion production will move the system towards chemical 
equilibrium by lowering the (N/Z) of the daughter system constrained by total charge con- 
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Fig. 25. Upper left panel: Excitation function of the 4-7r-integrated ratio of ir~ /ir + yields in 
central Au+Au collisions. The experimental data are joined by a least squares fit of the function 
Co + C-\(E / A)~ l excluding the lowest energy point. The IQMD SM prediction (triangles) is also 
given. Upper right and lower left panels: the N/Z dependence at 1.5^4, respectively OAA GeV of 
the 7r~/-7r + ratio. The solid lines are least squares fits of linear or quadratic (N/Z) dependence. 
Lower right panel: same as lower left panel, but for filtered data. 



servation [72]. However, the linear behaviour of IQMD at OAA GeV, where pion emission 
is a modest perturbation, is less trivial, all the more as it disagrees with the non-linear 
behaviour of the data. 

A systematics of tt~ /tt + ratios was first established for inclusive reactions [51] at 0.8A GeV 
beam energy using various, also asymmetric systems. In Fig. 26 we reproduce these older 
data which were plotted as a function of an estimated [51] 'fireball' (N/Z) composition. 
Due to the limited accuracy, both linear and quadratic (N/Z) dependences are compatible 
with these inclusive data. Our data point for Au+Au at the same energy, but for a central 
collision selection, is also shown in the figure and is perfectly compatible with the linear 
extrapolation. 
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Fig. 26. Left panel: 7r~ /ir + ratio versus N/Z of the 'fireball' measured [51] in various inclusive 
heavy ion reactions at 0.8^4 GeV (full circles). The solid (dashed) curve is a linear (quadratic) 
least squares fit to the data constrained to be one at N/Z = 1. The data point (open circle) 
from the present work holds for Au+Au and was not included in the fit. Right panel: n~ /ir + 
ratios versus beam energy obtained in transport calculations [70] for the system 132 Sn+ 124 Sn 
(N/Z = 1.56), using two options for the symmetry energy, Esy-b, dashed, and Esy-a. The ratios 
obtained from the present Au+Au data (solid circles) by linear extrapolation (from N/Z = 1.494) 
are shown for comparison. 

In ref. [70] Bao An Li has performed more advanced calculations for 132 Sn+ 124 Sn at OAA 
GeV using two different options for the asymmetry energy, Esy-a, a version that increases 
linearly with the density, and Esy-b, an exotic variant that after an increase at lower 
density bends back down to cross zero again at p = 3po- The changes in the n~ /n + ratio 
that these very different alternatives induce are shown in Fig. 26 which reproduces results 
from Fig. 2 of ref. [70]. 

The N/Z ratio of the rare isotope beam combination, N/Z = 1.56 is not very different 
from that of 208 Pb+ 208 Pb, N/Z = 1.536 (which would be a more readily available bigger 
system). For the reaction 197 Au+ 197 Au studied here, N/Z = 1.494, we anticipate, by linear 
extrapolation to N/Z = 1.56, ir~ /ir + ratios that have been added to the figure. With this 
addition the problematics already known from the total pion multiplicity studies show 
up. First, the difference predicted from the calculation [70] between two rather extreme 
options is on the 10% level and hence on the order of the present experimental accuracy, 
and, second, none of the two predictions follows the data. Similar conclusions follow if we 
take the calculations in ref. [71]. 



9 Azimuthal correlations ('flow') of pions 



Owing to collective flow phenomena, discovered experimentally in 1984 [73,74], it is possi- 
ble to reconstruct the reaction plane event-by-event and hence to study azimuthal correla- 
tions relative to that plane. We have used the transverse momentum method [75] including 
all particles identified outside the midrapidity interval \y \ < 0.3 and excluding identified 
pions to avoid autocorrelation effects. 
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We use the well established parameterization 

U = (7, fa) ; u t = fitl 



dN 



vq[1 + 2vi cos(0) + 2v 2 cos(20)] 



u t du t dyd(f) 

v = v (y,ut) ; «i = v^y^ut) ; v 2 = v 2 (y,u t ) 
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=< cos(20) > 



where is the azimuth with respect to the reaction plane and where angle brackets indicate 
averaging over events (of a specific class). The Fourier expansion is truncated , so that 
only three parameters, vq,v\ and v 2 , are used to describe the 'third dimension' for fixed 
intervals of rapidity and transverse momentum. Due to finite-number fluctuations the 
apparent reaction plane determined experimentally does not coincide event-wise with the 
true reaction plane, causing an underestimation of the deduced coefficients V\ and v 2 which, 
however, can be corrected by studying sub-events: we have used the method of Ollitrault 
[76] to achieve this. The finite resolution of the azimuth determination is also the prime 
reason why the measured higher Fourier components turn out to be rather small. 




Fig. 27. The four quadrants in a plane, xy, 
transverse to the beam direction, z. The flow 
direction, x, is indicated. 



Alternatively to the three Fourier coefficients, one can introduce the yields Qi, Q 2 , Q3, Q4 
in the four azimuthal quadrants, see Fig. 27, of which only three are independent (on the 
average over many events) due to symmetry requirements 



Q2 — Qi 

Q = Qi+Q2 + Q 3 + Qa 

Q24 — Q2 + Qi 



The two equivalent triplets 

u , v u v 2 < — ► Qq, Qi, Q 24 

are related by 
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These relations show that v± is a dipole, while v 2 is a quadrupole strength. Statistical 
(count rate) errors can be deduced with use of elementary algebra. The quadrant formula- 
tion also has some advantage in assessing small apparatus distortions leading to systematic 
errors, as we shall see. The fact that v±, as well as v 2 , are found to be non-zero, is generally 
called 'flow' in the literature and in particular the first Fourier coefficient is taken to be a 
measure of 'directed flow', while the second Fourier coefficient has been dubbed 'elliptic 
flow'. 

Recently more general methods based on Lee- Yang zeros have been proposed [77] and 
used [78] to isolate 'true' collective flow from other effects that might influence azimuthal 
anisotropies. An analysis of some of our pion data with this method will be published 
elsewhere [78]. Here we take the point of view that the most likely 'non-flow' correlation, 
the A-correlation between nucleons and pions is part of the physics relating the pion flow 
to the nucleon flow, a physics that is implemented in the IQMD code that we use to try to 
understand the 'flow' observables. Removing this correlation from the data would require 
removing it from the simulation as well when comparing. In view of the high statistics 
required to do this with sufficient accuracy we prefer in this survey to stick to the simpler 
'standard' method. This allows us also to compare with older data where possible. 

Azimuthal anisotropies of pion emission in heavy ion reactions were first reported in the 
refereed literature by Gosset et al. [79] using the DIOGENE setup in the reaction Ne+Pb 
at 0.8A GeV. The authors attributed their observations to a target (Pb) shadowing effect. 
This study concerned primarily what is now termed the v\ component. This V\ compo- 
nent was studied in more detail in ref. [80] for the reaction of Au+Au at 1.15^4 GeV. 
The fact that in sufficiently non-central geometries the n + mesons had the opposite di- 
rected flow than the protons (termed 'antiflow') was reported supporting earlier theoretical 
predictions [81,33]. It was also observed that tt~ mesons had a different flow and it was 
concluded [80] that 'the differences between the behavior of the tt + and the tt~ suggest that 
further consideration of this phenomenon is needed'. 

The observation of a v 2 component was reported by the KaoS [82] and the TAPS [83] 
collaborations in 1993. A first shot at a v 2 systematics (Bi+Bi at 0.4, 0.7 and 1.5 A GeV) 
was made in ref. [84]. A qualitative interpretation of the anisotropies making use of the 
transient vicinity of spectator matter close to the participant fireball as a 'clock' was 
presented in ref. [85]. In terms of the azimuthal quadrants just introduced, the authors 
considered the ratios Qi/Q 2 , Q3/Q2 an d Q1/Q3 which mix the two Fourier components 
V\ and v 2 . 

A direct quantitative comparison of our present data and simulations with all these earlier 
observations is not straight-forward, as we shall see. 

For completeness we mention here that pion flow has also been studied at much higher 



beam energies, 1L4 GeV (AGS [86]) and 40A GeV, as well as 158A GeV (SPS [87]). These 
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incident energies are rather far apart, making it difficult to assess at present the gradual 
evolution of the pion flow observable and the associated physics changes. 

As pion 'flow' turns out to be small, when compared to nucleonic flow, its observation 
requires a high degree of systematic and statistical accuracy. Presently the latter is diffi- 
cult to achieve in transport code simulations if one has the ambition to predict vi(y,u t ) 
and V2{y,u t ) in their full two-dimensional glory. Despite our computational efforts the 
simulations are still plagued with statistical errors that are a factor 3-5 larger than the 
experimental ones. Therefore, in the sequel we shall use several flow characterizers, < v\ >, 
< viUto >=< u x o >, < t>2 > and < —V2uf >=< Uy X >=< Uy Q — ul > averaged (<>) over 
more or less large regions of phase space (besides the averaging over events of a specific 
centrality class). Although we present dipole and quadrupole flow in two separate subsec- 
tions, they should be considered as two sides of the same (rather complex) phenomenon. 
In the sequel averaging brackets will be omitted, but will be implied. 



9. 1 Directed flow v\ 



First, a few technical points will be mentioned. The near-equivalence of using the 'cos(0)' 
method, eq. , or alternatively the 'quadrant' method, as well as the effect of the resolution 
correction are shown in Fig. 28. 
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Fig. 28. Rapidity dependence of the directed ir~ flow in the reaction Ru + Ru at 1.5^4 GeV. 
The centrality is 0.25 < &o < 0.45. Left panel: Comparison of the 'cos(i^)' and the 'quadrants' 
methods. Right panel: Comparison of the resolution corrected data with the uncorrected data. 

Following [76] the correction factor is given by 

v n = v'J < cos(nA$) > 

where v' n is the (uncorrected) measured value and A$ = — $^ is the azimuthal angle 
between the true and the measured reaction planes, and respectively. Due to the 
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Table 1 



Reaction plane resolution correction factors for E/A = 1A GeV. 



< bo > 


Au+Au 


Ru+Ru 


Ca+Ca 


Au+Au 


Ru+Ru 


Ca+Ca 






cos(A<I>) 






cos(2A$) 




0.100 


0.852 


0.707 


0.561 


0.583 


0.365 


0.216 


0.167 


0.910 


0.786 


0.609 


0.708 


0.471 


0.259 


0.360 


0.963 


0.892 


0.698 


0.860 


0.663 


0.353 


0.502 


0.958 


0.894 


0.699 


0.845 


0.670 


0.354 



large acceptance of FOPI the correction factors are essentially given by the 'natural' finite 
number fluctuations rather than by apparatus limitations. Typical values for different 
average centralities at an incident beam energy of 1A GeV are shown in Table 1. The 
effect of the correction is smallest when nucleonic flow is large, i.e. in the third centrality 
bin of the heaviest system. Small inhomogeneities in the laboratory azimuthal acceptances 
were found to have negligible influence on the flow results. 



For the symmetric systems that we study here, V\(y) should be asymmetric with respect 
to midrapidity (yo = 0), in particular i>i(0) should be zero. As can be seen in Fig. 28, v±(y) 
does not cross the origin of the axes. A systematic study of this mid-rapidity offset showed 
that it depended on particle type, centrality and system size in a way suggesting that it 
was correlated with the track density difference in the 'flow' quadrant Qi and the 'antiflow' 
quadrant Q%. While this could be simulated using our GEANT based implementation of 
the apparatus response, a sufficiently accurate quantitative reproduction of the offset at 
mid-rapidity was not achieved. We therefore opted for an empirical method to correct 
for the distortion, using the very sensitive requirement of antisymmetry with respect to 
midrapidity. 

One can show that a good first order correction simply consists in shifting Vi(y) down 
by a rapidity independent correction Av± until it crosses zero exactly at mid-rapidity. 
Essentially V\ ~ (Qi — Qz)IQo assume we correct Qi by replacing it by Qi + cQi where 
c is a constant not depending on transverse momentum or rapidity (some global loss in 
the flow quadrant). The correction to V\ then is Avi ~ cQi/Q ~ c/4 where we use the 
approximation Qi ~ l/4Qo- After trying different Ansatzes for the correction c = c(y,u t ) 
which could be in principle a 2-dimensional function, we ended up using the rapidity 
independent Ansatz c = v o c + V\ c * u t which grows linearly with ut (or pt); Vq c andt>i c are 
parameters. The track density distortions are assumed to grow linearly with the density. 
The quadrant Q% is left unchanged (only relative corrections matter), Qi is replaced by 
Qi (1 + c) and Q24 by Q24 * (1 + 0.55c) implying that in first approximation it is sufficient to 
use Q 2 ~ 0.55*(Q! + Q 3 ). Using the factor 0.55, (±0.05), in the Q 2 4 correction, rather than 
0.50, the value for isotropic emission, allows roughly for globally enhanced out-of-plane 
emission in the present energy regime. The values of vq c and v\ c are fixed from the two 
mid-rapidity conditions f 1(2/0 = 0) = and u x (y = 0) = 0. This is done separately for 
each particle, centrality and beam energy. The systematic error of v\ after the distortion 
correction is assessed to be 0.007 or 7% (whichever is larger) from the uncertainty of 
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Fig. 29. Rapidity dependence of the directed flow, i>i(yo) ; for Au+Au collisions at \.hA GeV 
with centralities b < 0.25, 0.25 < b < 0.45, 0.45 < b < 0.55 (from left to right). The data are 
taken in the interval 1.0 < uto < 4.2. (Red) squares: ir + , (blue) diamonds: ir~ . The solid curves 
are least-squares fits, see text. 

the mid-rapidity offset determination. The systematic error of v 2 (and Uy X )) is (7 — 10%) 
and was estimated by varying the coefficient for the Q24 correction between 0.5 and 0.6. 
These systematic errors should be kept in mind when inspecting the figures presenting 
our flow data, as these contain only the statistical errors. In comparative cases, such as 
assessing isospin differences between the systems Zr+Zr and Ru+Ru the systematic errors 
are expected to cancel to some degree. 

We are now in a position to review a representative sample of the rich flow data that 
we obtained in the present work. We start with the well known 'S-shaped' curves for the 
rapidity dependence of flow, Vi(y ). Fig. 29 shows the data for Au+Au at 1.5^4 GeV and 
for three of our 'standard' centrality intervals. Several remarks can be made. Only the data 
for 7/0 < were actually measured, (anti) symmetry was used to infer the yo > behaviour. 
As bo is increased, the diagrams 'rotate' clockwise, the 7i + data always 'preceding' the n~ 
data. In the interval around bo = 0.5 this 'rotation' has moved into a new 'quadrant', the 
antiflow side, for the 7r + . It is clear from this observation that pion flow cannot be simply 
derived from a 'parent' A baryon flow assuming the latter to be equal to that of single 
protons, as was attempted in ref . [84] (for the elliptic part of the flow) . A key question of 
theoretical analysis will be whether the tt + /tt~ difference ('isospin differential flow') is just 
a Coulomb effect or rather necessitates a (nuclear) isospin effect in order to be reproduced. 

In Fig. 30 we present the transverse momentum dependence, Vi(u t o), for data integrated 
over the backward hemisphere. The switch of sign for 7r + near 60 — 0.5 is again visible. 
There is a puzzling weak wavy aspect, especially in the 7r + data that is at the limit 
of apparatus distortions that we cannot completely exclude. However, a simulation with 
IQMD, not shown here, predicts qualitatively similar, even more pronounced, structures. 

Although the detailed shapes of the flow data are seen to be complex, it is of some inter- 
est to try to characterize flow with just one parameter, thus easing the task of obtaining 
a survey of many system-energy-centralities. One such parameter, the midrapidity slope 
dv\(yo) I dyo\yo=o, stresses more the midrapidity region, while the alternative parameter, 
which we term 'large acceptance' flow is v\ (or u x q) averaged over a large region of mo- 



38 



Au+Au 1 .5A GeV 



- 1 1 ' 1 ' 1 ' - 
b < 0.25 


- 1 1 ' 1 ' 1 ' - 
0.25 < b < 0.45 


- 1 1 1 1 1 - 

^ n ,-, zcPhPhPiSt 

l-^mu 


□ nn D n 

□□□□□□q „ 

°°oooooc.o* yD Qg| 

- , 1,1,1,- 


□ □□□ 

D n nnn nnn D 

5 

- , 1,1,1,- 


o o J-ITS 
Ooo° i 

0.45 < b < 0.55 



J i I i I i I I I i I i I i JZ i I i L_ 

1 2 3 1 2 3 1 2 

u to u to u to 



Fig. 30. Transverse momentum dependence of the directed flow, vi(uto), for Au+Au collisions 
at 1.54 GeV with centralities b < 0.25, 0.25 < b < 0.45, 0.45 < b < 0.55 (from left to right). 
The data are taken in the rapidity bin —1.8 < yo < 0.0. (Red) squares: ir + , (blue) diamonds: ir~ . 




Fig. 31. Centrality dependence of the midrapidity slope dvi/dyo at various indicated beam 
energies in the reaction Au+Au. Open squares joined by dashed lines: 7r + , full diamonds: ir~ . 

mentum space and stresses more the regions closer to target (or projectile) rapidity since 
V\ is zero at y = 0. 

A survey of midrapidity slopes in Au+Au reactions as a function of centrality is shown for 
various indicated beam energies in Fig. 31. The slopes were determined by least squares 
fitting with the polynomial t>i + Vuy + v^y^ in the rapidity range —1.8 < y < (the 
data are averaged over the interval 1 < u to < 4.2). The quality of the fits can be visualized 
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Fig. 32. Centrality dependence of directed flow, v±, for various indicated system masses 
and incident beam energy 1.5^4 GeV. The flow is averaged over ir + and tt~ . The values for 
A p + A t = 96 + 96 are an average between Zr+Zr and Ru+Ru. 

in Fig. 29 (smooth curves). The constant t> 10 , which should be zero, accounts for remaining 
uncertainties of the offset correction described earlier, i> 13 is necessary because, obviously, 
vi(yo) is not linear over the extended rapidity range and, finally, vn is identical to the 
slope dvi(yo) / dyo at yo = 0. The 7r + slope data are seen to vary moderately with energy, 
the gap to the n~ slopes widens at the lower energies. 

System size dependences and the system isospin dependences are shown in terms of the 
large acceptance flow in Figs. 32 and 33. For the averaging we choose the intervals — 1.8 < 
y < and 0.8 < u t0 < 4.2 which are well covered by our setup for all measured system- 
energies. Although the data were obtained in the backward hemisphere we give the values 
in the forward hemisphere (which are opposite in sign) so that positive (negative) values 
mean 'flow' ('antiflow') analogue to the midrapidity slopes. For the size dependence study 
we have removed the isospin difference by averaging over ir~ and 7r + data and also, for 
mass 96 + 96, by averaging over the two systems 96 Zr + 96 Zr and 96 Ru + 96 Ru. The system 
size dependence of directed flow is seen to be complex, in contrast to elliptic flow (see next 
subsection) . 

The last two mentioned systems were also the basis for an isospin dependence study, see 
Fig. 33. We find for v\, as well as for u x q (both in large acceptance), an indication that the 
7r~ — 7r + difference is slightly larger for Ru+Ru than for Zr+Zr. The Coulomb potential 
difference between the two systems is expected to be close to 10%, a difference that does 
not seem to account quantitatively for the observations. Since N/Z is larger for the Zr+Zr 
system, this can only be understood if there is a nuclear isospin effect opposite in sign to 
the Coulomb effect. This seems to be also suggested by the theoretical investigations of 
Qingfeng Li et al. [44]. We shall come back to this when adding later also information on 
isospin differential elliptic flow (Fig. 40) which is easier to parameterize. 

Before trying to assess the sensitivities of pion flow data to theoretical input, we shall 
compare our data with those of ref. [80] which were taken at 1.15A GeV for the system 
Au+Au. Such comparisons are not trivial, as the system-energy, the centrality, the chosen 
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Fig. 33. Comparison of the centrality dependence of directed flow {v\ bottom and u x q top) in 
the systems 96 Zr+ 96 Zr (left) and 96 Ru+ 96 Ru (right). Open squares joined by dashed lines: 7r + , 
full diamonds: 7r~. 

specific flow-describing observable, and, last but not least, the phase space covered, must 
be 'aligned'. Fig. 34 presents in three panels the data relevant for this comparison in 
terms of the centrality (6 ) dependence of the mid-rapidity slope du x o/dyo. As we have 
no measurement at 1.15^4 GeV, we show in two of the panels our data at 1.0A and 1.2 A 
GeV. The analysis of the 1A GeV data has been extended to higher 6 to allow for a more 
complete comparison, but one can say that there is no dramatic difference between the 
two energies. The Kintner data have been converted to the same, scaled, axes using the 
information from Fig. 3 of ref. [80]. Similar to the slopes of vi(yo), our mid-rapidity slopes 
were determined from a least squares fit of u x o(yo) to the polynomial Uioo + ^xoil/o + ^x03l/o 
in the backward hemisphere —1.8 < yo < 0. At midrapidity du x o/dyo = u x qi, while the 
constants w^oo an d u X Q3 take care of the uncertainty of the offset correction and of non- 
linearities in the rapidity dependence, respectively. The quality of these fits is generally 
excellent and similar to those shown in Fig. 29 for t>i(|/o)- 

Interpolating our data at 1A GeV, we find compatibility with the data of Kintner et al. at 
b = 0.6, but for lower 6 the evolution is different. The crossing point to 'antiflow' of the 
earlier data is seen to occur roughly around b = 0.2 for n + , whereas our data suggest this 
to happen for 6 — 0.45 — 0.55. This difference exceeds the uncertainty of the centrality 
determination in both experiments. For 7r~ the data of ref. [80] are very close to the no- 
flow axis (given the error bars) preventing a reliable determination of the crossing point, 
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Fig. 34. Centrality dependence of the midrapidity slope du x o/dyo in the reaction Au+Au 
(squares and dashed lines: tt + , full diamonds: ir~). The left panels represent the data of the 
present work ('FOPI') for 1.CL4 and 1.2A GeV beam energy with polynomial fits and a cut 
1.0 < uto < 4.2. The right panel is adapted from ref. [80] for 1.15^4 GeV. 

which appears in our data to be shifted by about 0.1 units relative to the ir + curve. Also, 
our data show more 7r~ flow in the maximum near bo = 0.25. The midrapidity observable 
du x o/dyo depends on the transverse momentum range covered: in our case sharp cuts were 
applied (1.0 < u t0 < 4.2, or 0.10 < p t < 0.43 GeV/c at E/A = 1A GeV). The low p t 
limitations of the data were not discussed in ref. [80] . 



Moving now to a comparison with simulations, we have chosen the large acceptance flow, 
1 2/0 1 < 1-8 and 0.8 < u t o < 4.2, for this purpose as it requires less events for a given 
statistical accuracy. The comparison is presented in Fig. 35 for the reaction Au + Au at 
1.5 A GeV where the experimental data for Vi(bo) and u x o(bo) and the isospin differential 
flow Du x0 = u x o(tt~) — u x o(n + ) are framed by calculations with a soft EOS (left panel) and 
a stiff EOS (right panels). The stiff EOS predicts statistically significantly higher values of 
both Vi and u x q and a transition to antiflow in more peripheral collisions. This sensitivity 
was noted earlier [9]. In the 'standard' version of IQMD we use, the stiff EOS is closer to 



42 



IQMDSM FOPI IQMDHM 



0.05 

o 

x 
=J 

Q 

0.00 

0.05 
=3° 0.00 
-0.05 



T 

i 

I 
i 


1 

Au+Au 1 .5A GeV 


I T - 
I A 

5 








— i — i — i — — i — i — i — — i — i — i — [ 


— i — i — i — — i — i — i — — i — i — i — \ 


• — i — i — i — — i — i — i — — i — i — i — \ 


: <\ = 


— ^a. — 
" s ^\ 


S- -n 

a- -a. ^» 


_ ^ N 

\ _ 
\ 

\ 

1 1 1 1 1 1 h^H \ 


; h : 


■ x ^ \ 

\ 

\ 

^ I 




— i — i — i — — i — i — i — — i — i — i — [ 


— 1 — 1 — 1 — — 1 — 1 — 1 — — 1 — 1 — 1 — [ 




\ 




- ^ 

\ 

\ 

\ 

^ - 


\ 

□ 

^ ♦ 77~ - 

□ - - 77 + 
i i I i i i I i i i 


X - 

\ 

\ 

\ - 

^ - 
1 , , , 1 , i , 



0.2 0.4 0.6 0.2 0.4 0.6 0.2 0.4 0.6 

b b b 



Fig. 35. Centrality dependence of large acceptance (see text) directed flow in the reaction 
Au+Au at 1.5 A GeV (squares and dashed lines: ir + , full diamonds: 7r~). Going from left to right, 
the middle panels ('FOPI') represent the measured data for v%, u x o and the isospin differential 
flow Du x q. The left (right) hand panels are the simulated data using a soft (stiff) equation of 
state. The dotted horizontal lines are merely references to guide the eye. The result of a linear 
least squares fit to the experimental isospin-differential flow (top middle row) is repeated in the 
adjacent upper panels. The dotted curve in the lower left panel results from a calculation (for 
7r~) with the phase shift prescription. 

the data, especially for central (bo < 0.2) collisions. However, even the stiff EOS predicts 
a transition to antiflow at smaller bo than our data. This conclusion is not changed if we 
shift to the phase shift prescription for the A baryon lifetime although this option changed 
the pion multiplicities significantly (section 6.2). In the lower left panel of Fig. 35 we show 
also a calculation (dotted curve) with the phase shift prescription for 7r~. The statistical 
errors (not shown) are similar to those of the standard calculation. 

The systematic difference between ir~ and ir + flow is reproduced and seems to be unaffected 
by the EOS. Since we have not used an isospin dependent EOS, any isospin effect in the 
calculation, besides Coulomb fields, would have to result from the 'cascade' (i.e. collision) 
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Fig. 36. Rapidity dependence of the elliptic flow, V2(yo), for Au+Au collisions at 1.5A GeV 
with the indicated centralities. The data are taken in the interval 1.0 < uto < 4.2. The solid 
curves are least squares fits of the two-parameter function cq + C2]Jq- (Red) squares: 7r + , (blue) 
diamonds: 7r~. 

part of the model. Past experience has shown that cascade effects influence directed flow 
only moderately. The upper middle panel shows that Du x o(b ) varies linearly: the fit 
to the difference data is also plotted in the panels showing Du x $ from the calculations. 
The calculated Duq are similar to to the experimental values but, despite the statistical 
limitations, seem to indicate a small surplus at intermediate bo, possibly a consequence of 
the missing nuclear isospin mean field in the simulation. The effects are small and hence 
difficult to assess quantitatively in a convincing way. 

9.2 Elliptic flow V2 

The presentation of our elliptic flow data follows in many ways the scheme of the previous 
subsection on directed flow. We start with aspects of rapidity and transverse momentum 
differential flow in the reaction Au+Au at 1.5 A GeV, then switch to a systematics of 
beam energy, system size and system isospin dependences. Then we compare our data to 
earlier measurements where possible and finally compare a subset of the data to theoretical 
simulations. 

In Fig. 36 we show the rapidity dependence of elliptic flow, ^2(3/0)5 f° r Au+Au collisions 
at 1.5 A GeV for various centralities. Again, the flow of ir + is different from that of the ir~, 
the difference increasing with impact parameter. The rapidity dependence is relatively flat: 
elliptic pion flow is not limited to a narrow mid-rapidity interval. There seem to be some 
weak structures (we remind however that only the backward hemisphere was covered by our 
analysis). These marginal structures can be evened out by the two-parameter polynomial 
fits (see the caption) on a level that is close to the statistical freedom and therefore will not 
be discussed further. The transverse momentum range is again sharply cut in the limits 
1 < uto < 4.2. Note that we always plot (— v?), rather than (+V2). The sign of v% tells us 
that there is a surplus of emitted pions in the direction (y) perpendicular to the reaction 
plane (zx). 

The transverse momentum dependence of elliptic flow, v 2 (u t o), integrated over the back- 
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Fig. 37. Transverse momentum dependence of the elliptic flow, ^2(^*0); f° r Au+Au collisions 
with indicated centralities. The rapidity interval is —1.8 < 7/0 < 0. The solid curves are least 
squares fits of the polynomial ciu t0 + C2U^ + c^uf to guide the eye. (Red) squares: ir + , (blue) 
diamonds: 7r~. 
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Fig. 38. Centrality dependence of large acceptance elliptic flow vi (—1.8 < yo < and 
1 < uto < 4.2) at various indicated beam energies in the reaction Au+Au. Open squares joined 
by dashed lines: tt + , full diamonds: tt~ . 

ward hemisphere, is shown in Fig. 37. It rises with transverse momentum but seems to 
saturate beyond u t o = 3 (or p t = 0.375 GeV/c at 1.5 A GeV). The difference between n + 
and 7r~ flow is small but statistically significant for intermediate centralities and momenta. 

Beam energy dependences, system size and system isospin dependences are established 
again in terms of large acceptance flow. The outcome is summarized in Figs. 38, 39, 40. 
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Fig. 39. Centrality dependence of elliptic flow, V2, (—1.8 < yo < and 0.8 < u t o < 4.2) 
for various indicated system masses. The flow is averaged over tt + and ir~ . The values for 
A p + At = 96 + 96 are an average between Zr+Zr and Ru+Ru. The straight lines are linear 
least square fits constrained to vi = for bo = 0. 



Like vi, t>2 evolves only slowly with energy (Fig. 38), it + elliptic flow always being somewhat 
larger then 7r~ flow. The size dependence of v<i is rather pronounced (Fig. 39) and regular 
in contrast to v\ (Fig. 32). Comparing the systems 96 Ru+ 96 Ru and 96 Zr+ 96 Zr, Fig. 40, one 
can make the same comments as for v%, Fig. 33. The 7r + — 7r~ difference is larger for the 
system with the smaller N/Z. If we parameterize the bo dependence by (one-parameter) 
straight lines constrained to pass the origin, we see already by eye inspection (see Fig. 40) 
that the difference between the two systems in terms of isospin differential flow is larger 
than expected naively on account of the 10% Coulomb potential difference. Using the 
slopes of these least squares fitted lines, we can put this in a quantitative form. For u 2 the 
slope difference is 0.0127±0.0033 (Zr+Zr) versus 0.0221 ±0.0029 (Ru+Ru), and 

for u 2 y - u 2 x it is 0.054 ± 0.020 (Zr+Zr) versus 0.106 ± 0.010 (Ru+Ru). These observations 
suggest that a more extensive theoretical analysis of such data could help establishing 
constraints on isospin dependences of high density mean fields. 

We close this subsection by comparing our data first to KaoS data [84] and then to IQMD 
simulations. KaoS published elliptic flow data for Bi+Bi at 0.4, 0.7 and 1A GeV. As we 
do not have data at 0.7 A GeV, and since the statistical significance of our measurements 
at 0AA GeV is rather modest, we limit ourselves to a comparison at 1A GeV ignoring the 
difference between the systems Bi+Bi and Au+Au. The comparison is shown in Fig. 41 
for three centralities which were dubbed MUL2 (6 = 0.66), MUL3 (6 = 0.49) and MUL4 
(bo = 0.29) in ref. [84]. (The centrality selection of KaoS is based on charged particle 
multiplicities registered in a 96- units hodoscope at polar angles between 12° and 48°.) 
Again, we try to align the experimental conditions as much as possible, choosing similar 
nominal centralities and rapidity cuts for our analysis. The plotted KaoS data were inferred 
from Tables 1 and 2 of ref. [84] which tabulate P 2 — 2f 2 values uncorrected for resolution 
and then corrected using correction factors from Table 4 of ref. [88] . The average bo were 
inferred from Table 2 [88] using b = b/R with R = 13.468 fm for Bi+Bi. The shown KaoS 
data can be smoothened by a polynomial fit (see caption Fig. 41). 
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Fig. 40. Comparison of the centrality dependence of elliptic flow (v% bottom and Uy — u x top) 
in the systems 96 Zr+ 96 Zr (left) and 96 Ru+ 96 Ru (-1.8 < y < and 0.8 < u t0 < 4.2). Open 
squares joined by dashed lines: tt + , full diamonds: ir~ . The straight lines are linear least square 
fits constrained to v% = for bo = 0. 



We find an excellent agreement of the data for 7r~ flow, but not for tt + flow at b = 0.29 and 
0.49. This is somewhat puzzling as the authors of ref. [84] stress that their data show no 
significant difference in the azimuthal emission pattern of positively and negatively charged 
pions. Under the conditions chosen to generate Fig. 41 (i.e. limiting the rapidity range to 
1 2/0 1 < 0.5) we indeed find that tt + and tt~ flow are very similar, although integration over 
a larger yo interval reveals a small enhancement of 7r + flow as we saw above. We also note 
that another statement in ref. [84] that 'the pion data indicate little dependence on the 
impact parameter' is not properly characterizing our data for bo < 0.5. 



In Fig. 42 we present a comparison of measured elliptic flow with results from simulations. 
Following a similar line as for directed flow (Fig. 35) we present large acceptance data 
for (—V2), (Uy — u 2 xQ ) and isospin differential {n + — n~) flow, D(iiy — u x0 ), in Au on Au 
collisions at 1.5 A GeV framed by simulations using a soft (left panels) and stiff (right 
panels) EOS. Qualitatively, the trends of the simulations are similar to those of the data, 
the stiff EOS again being somewhat closer to the data, although still a little on the short 
side. The systematic n + — n~ difference is well reproduced, see the upper row of panels 
in Fig. 42. As already noticed for directed flow (Fig. 35), the phase shift prescription for 
the A baryon lifetime leads to a moderate (and statistically marginal) modification of the 
predicted elliptic flow (dotted curve in the lower left panel for 7r~). 
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Fig. 41. Transverse momentum dependence of elliptic pion flow, V2(uto), for incident beam 
energies of 1A GeV and various indicated centralities b$. The FOPI data for Au+Au are compared 
to the KaoS data [84] for Bi+Bi. The smooth curves are least squares fits of the two-parameter 
polynomial ^22^0 + v 24 u to to the KaoS data guiding the eye. The data are taken in the rapidity 
interval |yo| < 0.5. 

10 Summary 



This work presents for the SIS energy regime of heavy ion collisions the most encompassing 
pion systematics available today. While some of the presented data are just a (needed) 
confirmation of earlier pioneering work, a significant number of observations are new. 

Among the data which urgently needed a confirmation, are pion multiplicities. One impor- 
tant conclusion from our present study is that two pion detection systems that were able to 
measure and publish pion multiplicities in heavy ion reactions in the 1A GeV regime with 
acceptances close to 4it, the BEVALAC Streamer Chamber and FOPI at SIS/Darmstadt, 
now offer highly consistent results (Fig. 15). In the larger framework of using heavy ion 
collision data to infer properties of nuclear matter far off the ground state, this should 
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Fig. 42. Centrality dependence of elliptic flow in the reaction Au+Au at 1.5A GeV (squares 
and dashed lines: 7r + , full diamonds: vr~). Going from left to right, the middle panels ('FOPI') 
represent the measured data for V2, Uy — v%, and the isospin differential flow D(v,y — u%)- The 
left (right) hand panels are the simulated data using a soft (stiff) equation of state. The dotted 
horizontal lines are references to guide the eye. The result of a linear least squares fit to the 
experimental isospin-differential flow (top middle flow) is repeated in the adjacent upper panels. 
The dotted curve in the lower left panel results from a calculation (for ir~) with the phase shift 
prescription. 

encourage future theoretical efforts to understand pion production in the 1A GeV regime 
on a level of 10%, the experimental uncertainty that is, presently still, typical for most 
absolute pion observables, and in particular the multiplicities. Extending the truly pio- 
neering Streamer Chamber data, our present data give information on two pion charges, 
7r + and tt~, and hence do not require non-trivial assumptions on the isospin dependences 
in order to deduce the multiplicity of all pions. 

Also, our systematics in terms of system sizes and system isospin is now enlarged. As a 
result of this improvement, we can now say that pion multiplicities per participant are 
not strictly a constant for a given incident energy, but show a measurable trend towards 
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smaller values as the participant number is increased, an effect that on the 10% level does 
not depend on how the system size is varied, be it by varying the overlap zone changing 
the centrality, or by modifying target and projectile (Fig. 18). 

Most of the progress since the eighties, however, rests on the detailed, centrality selected 
truly three-dimensional momentum space populations that have now become available. In 
contrast with just integrated 4n multiplicities which require large (and well understood) 
acceptances rather than a multitude of registered events, many subtle and relatively small 
effects, such as asymmetries in pion emission, required electronic devices such as FOPI 
to be able to handle the wealth of information connected with the variation of energy, 
centrality, system size and system isospin. Complemented with the studies of the TAPS 
and the KaoS Collaborations there now exists a rather complete set of informations on 
pion emission. 

In two (scaled) dimensions, longitudinal rapidity, y , and transverse four- velocity, u t0 , the 
pions peak at mid-rapidity populating, especially at low energy, a rather broad phase 
space extending over at least twice the original rapidity gap. In these scaled units the 
distributions become significantly more compact at the higher end of the studied energy 
range (Fig. 6). 

Rather than presenting the longitudinal (beam) direction in terms of a rapidity distri- 
bution on a linear scale, and then switch to a logarithmic ordinate scale to present the 
transverse direction in terms of a transverse mass or momentum distribution, as is com- 
monly done in the literature, we have compared the two orthogonal directions more directly 
introducing the (one-dimensional) transverse rapidity distributions. The scaled variances 
of the transverse and the longitudinal rapidity distributions decrease significantly with 
increasing beam energy (Fig. 7), the transverse variance always being the smaller of the 
two. We have dubbed the ratio of the transverse to the longitudinal variance 'stopping' 
and find that this observable decreases steadily with increasing incident energy, although 
more slowly than the individual variances, which in a naive thermal model (equilibrium, 
no flow, no decay distortions) would relate directly to (kinetic) temperatures and would 
have to be independent of the direction. 

This 'pion stopping' qualitatively follows the 'nucleon stopping' [48] , the latter decreasing 
however faster with energy, probably a reflection, on average, of a less violent collision his- 
tory, due in part to a less perfect separation of spectator type influences. However, we also 
found evidence for spectator influence on the pions, although created particles originate 
in this energy regime from the overlap zone. When studying the size dependence of pion 
stopping, we found a remarkable phenomenon of two stopping branches: while stopping in 
the most central collisions is almost independent of the system size, very much in contrast 
to nucleon stopping [48], and even seems to increase slightly as the size decreases, an effect 
that is highly non-trivial, we find that stopping decreases when significant spectator matter 
is present (i.e. in half-overlap, b = 0.5, collisions). We associate this observation (Fig. 8) 
tentatively to a reacceleration of pions penetrating and rescattering in the fast specta- 
tor (in the c.o.m. system). Since there was no significant indication of two corresponding 
branches in pion production, one is tempted to conclude that pions are rescattered, rather 
than absorbed by spectator matter. 

Isospin dependence of pion stopping was not observed on a detectable level when comparing 
96 Ru + 96 Ru with 96 Zr + 96 Zr. 
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The partial transparency that is suggested by the observation that the 'stopping' ob- 
servable is always less than one, is also evidenced by a marked polar anisotropy of pion 
emission (Fig. 9) which again shows the two branches depending on the presence or ab- 
sence of spectators (Fig. 11). Anisotropics imply that extrapolations of midrapidity data 
to 47r that assume isotropic emission underestimate pion yields typically by a factor 1.2 
to 1.4. When we apply this correction factor to TAPS or KaoS midrapidity data, we find 
fair agreement with the data from large acceptance devices, however some problem cases 
remain (Figs. 16, 18). 

Our transverse momentum spectra, that agree in shape with KaoS (Fig. 22), and also 
with the IQMD simulation (Fig. 23), are characterized, besides the variance (or apparent 
transverse temperature) by a marked difference between n + and 7r~ that is suggested 
by the simulation to originate almost exclusively from the Coulomb fields (Fig. 23). The 
latter also account for the main features of the systematics of ratios of average momenta 
(Fig. 24). 

Whether the symmetry part of the ground state EOS influences the yield ratios of n~ to 
7r + in neutron rich nuclei is not clear. The naive expectation from the first chance isobar 
model that this ratio should grow quadratically with N/Z is generally not fulfilled, as the 
ratio tends to vary linearly with N/Z due to a partial thermalization effect, except at 
the lowest energy (OAA GeV). Simulations seem to miss this low energy trend, however 
(Figs. 25, 26). 

This work also provides a significant extension of data on pion azimuthal correlations, 
especially for the directed flow which is of special interest due to its apparent sensitivity 
to the EOS (Fig. 35) [9]. 

The rapidity dependence ^i(yo) has the S-shape familiar from the early days of flow mea- 
surements [89] and is well described in terms of an odd polynomial in y including just a 
linear and a cubic term (Fig. 29). In contrast to the ultrarelativistic regime [87] where the 
flow of pions appears to be mostly opposite to proton flow ('antiflow'), we find in the 1A 
GeV regime a more complex behaviour. First, for high centralities, b < 0.4, the flow of the 
pions follows the nucleon flow although it is considerably smaller in terms of v±, roughly 
by the ratio of pion to nucleon mass. For larger impact parameters there is a switch to 
antiflow (Figs. 29 and 34). Second, the flow of positively charged pions differs from that of 
negatively charged pions, in particular the switch to antiflow occurs at smaller centrality 
(Fig. 34). 

Although these features confirm qualitatively what was found in an earlier study by the 
EOS Collaboration [80], the quantitative agreement is not satisfactory (Fig. 34): whether 
the difference is due to different detector acceptances is not clear. 

The transverse velocity dependence Vi(u t o) (integrated over a large rapidity interval) is 
relatively flat with weak structures (Fig. 30). 

Despite theses detailed features which show that pion flow is a complex phenomenon, 
it is useful for systematics purposes to try to characterize flow by just one parameter. 
One alternative is the midrapidity slope dvi/dy (y = 0), for which we have established a 
systematics varying the incident energy (Fig. 31): this parameterization puts more weight 
on the midrapidity region. Another alternative, stressing higher rapidities, consists of the 
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large-acceptance value of v±, for which we have established size (Fig. 32) and isospin 
(Fig. 33) dependences. The isospin differential flow of pions is a systematic effect which 
is largest for Au + Au at the lower incident energies (Fig. 31). It seems to be dominated 
by Coulomb effects since it is quantitatively reproduced by our simulation for Au+Au at 
1.5 A GeV (Fig. 35) which did not include isospin dependent mean fields except for the 
Coulomb field. However, surprisingly, the isospin differential pion flow was found to be 
significantly smaller for 96 Zr + 96 Zr (N/Z = 1.40) than for 96 Ru + 96 Ru (N/Z = 1.18), 
Fig. 33. Quantitatively, it seems necessary to introduce a nuclear isospin field of opposite 
sign to the Coulomb part of the field in order to explain this feature. The simulation at 
1.5A GeV (Fig. 35) showed that the EOS influences pion directed flow in a measurable 
way, favouring the stiffer EOS in the IQMD version that we used. 

In the 1A GeV energy range pion elliptic flow like nucleonic flow is negative, i.e. preferen- 
tially out-of plane . The presence of spectator matter in the expansion phase of the fireball 
is probably the most important 'geometrical' difference to the ultrarelativistic regime [87], 
our stopping studies indicating that pion rescattering in spectator matter plays a role, 
as was already suggested in ref. [33]. We find that f 2(1/0) is rather flat over a very large 
range of rapidities (Fig. 36). In contrast to directed flow, |f 2 j increases nearly linearly 
with decreasing centrality in the range covered by our data (b < 0.55). The transverse 
four-velocity dependence v 2 (u t0 ) (Fig. 37) is quadratic (7r~) or linear (ir + ) for smaller 
momenta, flattening out at higher momenta. We confirm the tabulated KaoS data [84] 
for 7r~ (Fig. 41), but not for ir + . While the dependence on incident energy is relatively 
weak, Fig. 38, the system-size dependence is rather strong for all centralities, Fig. 39. As 
for directed flow, the system 96 Ru + 96 Ru shows a larger difference between 7r + and 7r~ 
elliptic flow than the system 96 Zr + 96 Zr, so the same comments apply. In comparing with 
the IQMD simulations we find experimentally a somewhat larger elliptic flow, even when 
using the stiff EOS (Fig. 42). The isospin differential flow in the heavy system Au+Au is 
again well reproduced (Fig. 40). 

At this time, however we do not wish to draw firm conclusions on the issue of the stiffness 
of the EOS from the pion flow, as we expect other features of the simulation, namely the 
treatment of the A baryon propagation in the medium, to influence the predicted flow. 
Shorter effective A baryon lifetimes in the model are expected to decrease the efficiency 
of the assumed dominant pion absorption mechanism via NA — > NN and hence the 
pion yields at freeze out, as well as the influence of the mean field in the time between 
elementary collisions, and hence the observed flow. This is an interesting aspect of in- 
medium physics in its own right and deserves theoretical investigations with the aim of 
reaching a consensus on how this problem should be treated properly in transport codes. 

Much work still needs to be done, to coordinate the transport theoretical efforts with the 
aim of achieving code independent conclusions. Such efforts were started in ref. [27] and 
continued on a larger scale in ref. [64] which also contains further literature citations (and 
transport codes) for the interested reader 

The request for reproduction of finer details, such as the full 3D features of the momentum 
space population under exclusive conditions, subtle system size and isospin effects, is pre- 
mature at this time as long as global features, such as production, stopping and integrated 
flow are not under sufficient control. This is a task for future theoretical work beyond 
the scope of the present experimental work. From the richness of our observations it fol- 
lows that this represents a challenge and hopefully will eventually contribute to a better 
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understanding of nuclear medium properties, which cannot simply be inferred from inco- 
herently superimposing experimental information on so-called 'elementary' hadron-hadron 
reactions. 
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